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ABSTRACT 


▼ 


▼ 


This document is the final report of analyses and testing of stability 
augmentation systems accomplished under NASA- Langley Research Center Contract NAS1- 
10885 and is intended to be used as a working reference in future program activities. 
Section 2 describes the active flutter suppression analyses conducted on the NASA 
l/l7 scale supersonic transport wing model. Results predict a 16.7$ increase in 
model flutter true airspeed for Dr. Nissim's flutter suppression concept with out- 
board control surfaces and sensor locations. Mechanization of the flutter suppress- 
ion system on the wing model is discussed in Section 3. The analysis, design, fab- 
rication and testing of a subminiature electro-hydraulic actuator, weighing approx- 
imately 2 ounces and producing 4l in-lb. of torque, are described. Section k pre- 
sents the results of the ride control system synthesis for the NASA l/30 scale 
B-52 aeroelastic model. A ride control system design using elevator, flaperon and 
horizontal canards is described which reduces fuselage RMS vertical accelerations 
due to random gusts more than 30 $. 


RETRIEVAL REFERENCE WORDS: 
Aeroelastic Msdels 
Flutter Suppression Systems 
Ride Control Systems 
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1,0 INTRODUCTION 

This document is the final report of analyses and testing of stability 
augmentation systems accoirqplished from 25 May 1971 to 24 May 1972 under NASA-Langley 
Research Center Contract NASI- 10885 and is intended to be used as a working ref- 
erence in future program activities. 

Section 2 describes the flutter suppression system analysis conducted 
on the NASA l/l7 scale supersonic transport wing model, and Section 3 discusses 
the work that was accomplished on the mechanization of the flutter suppression 
system. The ride control system synthesis for the NASA B-52 aeroelastic model is 
discussed in Section 4. 

Each section is written to be independent of the other sections. The 
sections discuss not only the work accomplished, but also contain a general dis- 
cussion of the work remaining to con^lete the three work items. 
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2.0 SST WING MODEL FLUTTER SUPPRESSION SYSTEM ANALYSIS 

Evaluation of a flutter suppression system developed by Dr. Eliatau 
Nissim on the NASA one- seventeenth scale supersonic transport (SST) wing model 
is described in this section. This work was accomplished to demonstrate active 
flutter mode control on the model in the Langley transonic dynamics tunnel. 

2.1 Introduction 

An analytical study was conducted to determine performance of the 
flutter siqppression system on the SST wing model. A similar evaluation was con- 
ducted in 1970 with the flutter system on the 9^9-300 SST airplane configuration. 
Results of the airplane analysis are contained in Boeing document D 3 - 8 390-1 
(Reference l). This report is also contained in Section 2 of D3-8390-4 (Ref- 
erence 2 ) . 

Results of the airplane analysis showed a 28 percent increase in 
flutter speed at Mach 0.9 over the unaugmented airplane using midspan leading 
and trailing edge control surfaces. Based on these results, the midspan control 
surface location was chosen for the mechanization of the flutter suppression 
system on the model for the wind tunnel demonstration. The primary objective 
of the current analysis was to determine if the flutter speed could be increased 
sufficiently for wind tunnel demonstration of the flutter suppression system. 

The results show that a 12.3 percent increase in flutter speed at Mach 0.9 can 
be attained with these surfaces, with the feedback sensors located at the inboard 
edge of the wing strip. Although the flutter speed improvement is less than 
attained on the SST airplane, the midspan control surfaces will be used for the 
wind tunnel tests. 

The model analyses were conducted using equations of motion generated 
from generalized mass and stiffness data supplied by NASA. The Mach 0.9 equations 
were written to include the inboard and outboard, as well as the mid span, control 
surface locations to permit assessment of the system performance for all surface 
locations. The analyses described in Section 2.3 were conducted using the Mach 0.9 
equations, with the assunption that ideal actuators were driving the control sur- 
faces. To couplets the analysis, effects of non- ideal actuators must be determined 
and the system evaluated at Mach 0.6 and 1.2. 

2.2 Equations of Motion 

Equations of motion were developed for the wing model for Mach 0.6, 

0.9, and 1.2 using generalized mass and stiffness data supplied by NASA. The 
equations were written with wind tunnel velocity and fluid mass density as explicit 
functions to permit variations in dynamic pressure by varying either the velocity 
or mass density, or both. A 95 percent freon, 5 percent air environment was assumed 
for the wind tunnel fluid. 


The equations included six elastic degrees- of- freedom plus the two 
control surface displacements. The Mach 0.9 equations were subsequently revised 
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to include the inhoard and outboard leading and trailing edge control surfaces, 1 
with the actuator mass assumed at the inboard edge of each surface. Three sep- 
arate structural mass matrices were required to properly account for the actuator 
mass* 


In the equation generation, structural damping was assumed to be zero. 
Doublet- lattice unsteady lifting surface theory was used for the two subsonic con- 
ditions, and mach box theory for the Mach 1.2 condition, to obtain aerodynamic 
loading. The resulting coup lex matrices of unsteady aerodynamic coefficients 
were transformed through a curve fitting procedure to rational functions of the 
Laplace transform operator, S, with fourth order denominators. The equations 
were then rearranged to the form 

(s a [tA + pc,] + s£d+ fVCal + CK+ pv 2 c 3 "] 


where : 


[vQ, 

Cc,], 


pv 2 ZCd^Cs/cs + = {°} 


{ V s) } 
s 
p 

V 

W, [d] 

Ec 2 1, [C 3 ] 

C^Kjl 


= Elastic and control surface displacement degree s- 
of freedom 

= Laplace transform operator 
m Fluid mass density (95$ freon, 5$ air) 

= Velocity of fluid relative to the wing 
= Structural mass, stiffness, and damping 
= Aerodynamic parameters 
= Lift growth parameters 
= Aerodynamics parameters. 


Numerical values of the matrix elements for the three test conditions 
are presented in Section 2. 5* Locations of the control surfaces and doublet- 
lattice panels are shown in Figure 2.1. The sign convention used in the equations 
is: 


X - Positive aft 
Y - Positive outboard 
Z - Positive up 

Trailing edge surface displacement - Positive trailing edge down 
Leading edge surface displacement - Positive leading edge up. 
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The spanwise length of «n control surfaces included in the equations 
is 5*88 inches (11,76 percent of the wing semispan) with 20 percent chord width 
on all but the inboard and midspan leading edge surfaces. These two surfaces have 
a constant 3,65 inch width so they could be installed in the model without cutting 
into the aluminum alloy plate that forms the model elastic structure. 


2.3 Flutter Suppression System Evaluation 

The flutter suppression system analysis initially assessed the flutter 
speed improvement attainable with the control law identical to that used on the 
airplane analysis, A block diagram of the system is shown in Figure 2.2. Feed- 
back sensors were located at 30 percent and 70 percent of the wing chord at the 
center of the wing strip defined by the midspan leading and trailing edge control 
surfaces. The Mach 0.9 equations were expanded later to include the inboard and 
outboard control surface locations. This was done to permit an evaluation of 
system performance for the other control surfaces, including sensor location 
variations. 


All model analyses were conducted with the sixth elastic mode omitted 
from the equations of motion. This was done because of the unrealistic influence 
of this mode on the flutter characteristics. With the sixth mode included, the 
dynamic pressure at which flutter occurred for the unaugmented wing was about half 
the dynamic pressure at flutter with the sixth mode omitted. 

The flutter suppression system includes feedback variables whose ampli- 
tudes are proportional to displacements, but in phase with rates. These variables 
were generated for tjjie stability analysis using phase root locus to introduce the 
90° phase lead ) to displacements at all frequencies, without changing ampli- 

tude as a function of frequency. 

Dynamic pressure was varied by changing the fluid mass density in all 
the analyses except where the SAS performance was evaluated as a function of mass 
density. For this analysis, wind tunnel velocity was varied to vary dynamic 
pressure. 

2.3«1 Midspan Surfaces with Nominal Flutter Suppression System 

The flutter suppression system evaluation began with the midspan con- 
trol surfaces and the feedback sensors located at the center of the wing strip. 

The equations of motion used were for the Mach 0.9 condition without actuator mass 
included. 


A dynamic pressure root locus for the open and closed loop cases is 
shown in Figure 2.3. With the flutter suppression system loop open, the wing 
first elastic mode crosses the imaginary axis at about 136 lb/ft* dynamic pressure. 
Kie dynamic pressure at flutter with the loop closed is about 159 lb/ft^, an in- 
crease of 16.9 percent on dynamic pressure (8.1 percent in terms of wind tunnel 
fluid velocity). Figure 2.4 shows the open and closed loop damping ratios for the 
flutter mode as a function of dynamic pressure. 
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WHERE H, =• VERTICAL DISPLACEMENT AT 307 tt WING CHORD (POSITIVE UP) 
H 2 - VERTICAL DISPLACEMENT AT 70% WING CHORO (POSITIVE UP) 

OC = (POSITIVE LEADING EDGE UP) 

% = TRAILING EDGE SURFACE DISPLACEMENT 
S LE » LEADING EDGE SURFACE DISPLACEMENT 
c. = WING CHORD LENGTH AT SENSOR LOCATION 
b = C /a 


FIGURE 2.2 FLUTTER SUPPRESSION SYSTEM BLOCK DIAGRAM 








An analysis was conducted to evaluate system performance as a function 
of fluid mass density. Figure 2,5 shows the flutter mode damping ratio as a func- 
tion of true air stream velocity for variations in mass density. Effectiveness of 
the flutter suppression system as a function of fluid mass density is illustrated 
in Figure 2.6. The flutter speed improvement decreases with increasing mass 
density. 


2.3.2 Effects of Surface /Sens or Variations 

Since the flutter speed improvement attained on the model was consid- 
erably lower than expected from the previous airplane analysis, the Mach 0.9 equa- 
tions were modified to include the inboard and outboard control surfaces. The 
mass of each actuator was included to provide a more accurate representation of 
the model as it would be tested in the wind tunnel. 

The revised equations were used to assess the flutter suppression 
system performance attainable with each set of leading and trailing edge control 
surfaces, with the feedback sensors located at the inboard edge, center, and out- 
board edge of each wing strip as shown in Figure 2.7. The control law gains used 
in this analysis were the same as used previously for the midspan surfaces, center 
sensors configuration. 

Damping ratios of the flutter mode as a function of dynamic pressure 
for the three control surface configurations are shown in Figures 2.8, 2.9, and 
2.10. Each plot shows the effect of sensor location within the wing strip. Dyna- 
mic pressure root loci for the flutter mode are presented in Figures 2.11, 2.12, 
and 2.I3. 


Table 2.1 summarizes the results of this analysis. The outboard con- 
trol surfaces, with outboard edge sensor location, provides the most improvement 
in dynamic pressure at flutter, 36. 1 percent (or 16.7 percent increase in flutter 
speed assuming fluid mass density held constant). The midspan surfaces, inboard 
edge sensors were the next most effective surface/ sensor combination with 26.1 
percent increase in dynamic pressure (12,3 percent increase in flutter speed). 

TABLE 2.1 

EFFECT OF SURFACE/SENSOR LOCATION 
ON FLUTTER DYNAMIC PRESSURE 


SURFACE 

LOCATION 

DYNAMIC PRESSURE - psf 

OPEN 

LOOP 

CLOSED LOOP I 

INBOARD EDGE 
SENSORS 

CENTER LINE 
SENSORS 

OUTBOARD EDGE 
SENSORS 

INBOARD 

(BL 77.5 - 92.5) 

135 

159.5 

165. 

163. 

MID SPAN 

(BL 92.5 - 107.5) 

136 

171.5 

159.5 

147. 

OUTBOARD 

(BL 107.5 - 122.5) 

133 

lUl.5 

166. 

181. 
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FIGURE 2.5 VARIATION OF FLUTTER MODE DAMPING RATIO 
WITH TRUE VELOCITY (MID SPAN STRIP) 









FLUTTER 2.7 LOCATIONS OF SURFACES AND SENSORS 
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FIGURE 2.8 VARIATION OF FLUTTER MODE DAMPING RATIO 
WITH DYNAMIC PRESSURE (INBOARD SURFACES) 
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FIGURE 2.10 VARIATION OF FEXJTEER MODE DAMPING RATIO 

WITH DYNAMIC PRESSURE (OUTBOARD SURFACES) 
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FIGURE 2.11 DYNAMIC PRESSURE ROOT LOCUS 
(INBOARD SURFACES) 
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2.3.3 


SST Airplane Analysis 


A brief study was conducted to determine the effects of the flutter 
suppression system performance due to truncation of the mathematical model. The 
S T 969 - 3 OO configuration equations of motion were truncated to the first seven 
degree s-of- freedom from the 12 degrees-of- freedom used previously in the air- 
plane analysis. The control surface and sensor locations, and the control w, 
were identical to that used previously. 

The equations of motion were written with velocity and air density 
as explicit functions to permit varying forward velocity as a function of alti- 
tude^ constant Mach number to determine flutter speed. 

Figure 2.14 shows a dynamic pressure root locus for the two airplane 
flutter modes for the unaugmented full mathematical model and the seven degree- of- 
freedom equations. The mode damping ratios as a function of velocity (KCAS) are 
shown in Figure 2.15 with the flutter suppression system loop open and closed. 
Although thi basic airplane flutter speed was higher for the truncated mathematical 
model, the improvement attained with the flutter suppression system was about the 

same. 


2.4 


Remaining Work 


The flutter suppression system has been evaluated on the model Mach 
0.9 equations of motion with the assumption of ideal actuators driving the control 
surfaces. The effects of the actual actuator dynamics on system performance must 
be determined. To minimize mechanization difficulties, the midspan control sur- 
faces with inboard edge feedback sensors will be used during the wind tunnel tests. 

The flutter suppression system will also be evaluated at Mach 0.6 and 
1.2. The equations have been derived and once the Mach 0.9 analysis is conqpleted, 
the Mach 0.6 and 1.2 e valuations should be straightforward. 

When the control law to be mechanized on the model is finalized, sub- 
critical responses must be generated for correlation of wind tunnel test data with 
theoretical results. These responses should be generated on the most accurate 
mathematical model possible, including the non-ideal frequency measurir^ cir^H 
required by the control law mechanization as well as actual actuator dynamics an 
the effects of the actuators and hydraulic lines on the wing dynamic behavior, 
discussion of the control surface mechanization on the model is presented i 

Section 3. 


2.5 


Supporting Data 


Numerical values of the equation of motion coefficients are presented 
below. The Mach 0.9 equations contain three structural mass matrices correspond ng 
to the inboard, midspan, and outboard control surface lo^tions with 
at the inboard edge of each surface. The equations for Mach 0.6 and 1.2 contain 
coefficients for the midspan control surfaces only. 
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The motion variable vector, { > for the Mach 0*9 equations has 

dimension 12 x 1. The first six elements are the wing elastic modes; the seventh 
and eighth elements are inboard leading and trailing edge control surface displace- 
ments, respectively; the ninth and tenth elements are midspan leading and trailing 
edge control surface displacements, respectively; and the eleventh and twelfth ele- 
ments are outboard leading and trailing edge control surface displacements, res- 
pectively. For the Mach 0.6 and 1.2 equations, this vector is 8 x 1 with the first 
six elements the elastic modes and the seventh and eighth elements are the midspan 
leading and trailing edge control surface displacements, respectively. 

The structural mass, stiffness, and damping matrices are constant for 
the three wind tunnel conditions, except that the Mach 0.6 and 1.2 equations have 
aerodynamic coefficients for the midspan control surface location only. For these 
two test conditions, the first six and the ninth and tenth rows and columns of the 
mass, stiffness, and damping matrices are to be used. 
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O.OOOOE- 

01 

.0.00006- 

01 

O.OOOOE- 

01 

0.00006- 

01 

0 • OOOOE-0 1 

2. 5553 E 02 

■ - - : - 


— 








■ - — — 



NO. 

dr-8884 
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PAGE 
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_ MACH 0.9 

MATRIX *M* 12 x 12: MID SPAN SURFACE MECHANIZATION . 


ROW 1 


7.4782E 

02 

1.065 IE 

01 

-1.4161E 

01 

- 1 . 9078 £ 

01 

-3.9026E 

00 

2 • 99 30 E— 01 

3. 412 IE 

01 

-2.4087E 

02 

4.3332E 

01 

-1. 928 IE 

02 

4.1346E 

01 

-6 .2 867 EG 1 

ROW 2 

1.0651E 

01 

2. 7348E 

02 

-4.9 55 5E 

00 

—6. 1521E 

00' 

-3.6665E 

00 

2;0496'F~OrO"“ 

1. 1035E 
ROW 3 

01 

-3.9796E 

01 

2.9538E 

01 

- 1 .0223 E 

02 

3.6191E 

01 

—5. 0941 E 01 

-1.4161E 

01 

-4. 9555E 

00 

4. 8703E 

02 

7. 1 738 E 

00 

5.9995E 

01 

— 4.7 862 E 01 

-4.5802E 

01 

-1.8365E 

01 

-3.2175E 

01 

-6.427 IE 

01 

4.9542E 

00 

-5. 5438 E 01 

ROW 4 

"-U90 78E 

01 

-6. 1521E 

00 

7. 1738E 

00 

1.0278E 

02 ‘ 

-4.7197E 

00 

5.5 530fTO0“' 

-6.4176E 
ROW 5 

00 

1.3715E 

02 

-3. 1614E 

01 

8.9813E 

01 

-9.5785E 

00 

-2.6748E 01 

-3.9026E 

00 

— 3 • 666 5E 

00 

5.9995E 

01 

-4.7197E 

00 

3.5769E 

03 

-3.1 563 E 02 

-1.9838E 

02 

—9. 078 BE 

02 

- 1 . 8 63 2E 

02 

-6.0310E 

02 

-1.1190E 

02 

-1*041 5E 02 

ROW 6 

2.9929E- 

■01 

2.0495E 

00 

-4.7362b 

0 1 

5.5 53 OE 

00 

-3.1 563E 

02“ 

" 472452 E“0 2“ 

9.72C8E 
ROW 7 

01 

7.3279E 

02 

2.0682E 

02 

6.8965E 

02 

1.8020E 

02 

2.0033E 02 

3.4121E 

01 

1. 1035E 

01 

-4. 5802E 

01 

-6.4176E 

00 

-1.9838E 

02 

9.7 208 E 01 

3.4758E 

02 

O.OOOOE- 

■01 

O.OOOOE- 

01 

O.OOOOE- 

■01 

O.OOOOE- 

01 

O.OOOOE-Ol 

ROW 8 

— 2.4 087 E 

02 

-3.9796E 

01 

-1.B365E 

01 

1.3715E 

02 

-9.0788E 

02 

7.3279E 02 “ 

O.OOOOE- 
' ROW 9 

■01 

4.0587E 

03 

O.OOOOE- 

•01 

O.OOOOE- 

•01 

O.OOOOE- 

01 

O.OOOOE-Ol 

4. 3332 E 

01 

2.9538E 

01 

-3.2175E 

01 

-3. 1 614E 

Oi 

-1 .8632E 

02 

2.0682E 02 

O.OOOOE- 

01 

O.OOOOE- 

01 

2.8206E 

02 

O.OOOOE- 

01 

O.OOOOE- 

01 

O.OOOOE-Ol ' 

R0W10 

— 1.9281E 

02 

-1.0223E 

02 

-6.4271E 

01 

8.9813E 

01 

-6.0310E 

02 

6. 8965 E 02 

0. 00006- 
ROW 1 1 

01 

O.OOOOE- 

01 

O.OOOOE- 

01 

1.6917E 

03 

O.OOOOE- 

01 

O.OOOOE-Ol 

4. 1346E 

01 

3.6191E 

01 

4.9543E 

00 

-9.5785E 

00 

-1.1190F 

02 

1.8020E 02 

O.OOOOE- 

01 

0. OOOOfc- 

■01 

O.OOOOE- 

01 

O.OOOOE- 

01 

1.5570E 

02 

O.OOOOE-Ol 

ROW 12 

-6.2867E 

01 

-5.094 IE 

01 

-5. 5438E 

01 

-2.6748E 

01 

-1.0415E 

02 

2.0033E' 02 

O.OOOOE- 

01 

O.OOOOE- 

01 

O.OOOOE- 

01 

O.OOOOE- 

Ol 

O.DOOOE- 

01 

2. 55 53 E 02 




MACH 0.9 


MATRIX 5 M* 12 x 12: OUTBOARD SURFACE MECHANIZATION 


ROW 1 


7.86766 

0 2 

4.8375E 

0! 

6.3172E 

00 

- 3 * 0039 E 

01 -4.19436 

00 

2.Q902E 

01 

3.4121E 

01 

-2.40376 

02 

4*33326 

01 

-1.928 IE 

02 4.1346E 

01 

—6 .286/6 

01 

ROW 2 











4.8375E 

01 

3.0265E 

02 

3.7696E 

00 

-2.05266 

01 -6.07016 

00 

I • 8437 E"~ 

or 

1.1035E 

01 

—3 .97966 

01 

2.9538E 

01 

- 1 . 022 3 E 

02 3.6191E 

01 

-5.09416 

Oi 

ROW 3 











6.31726 

00 

3. 76966 

00 

4. 8063E 

02 

1.7884E 

00 5.23976 

01 

-6.8867E 

01 

-4.5802E 

01 

-1.8365E 

01 

-3.2175E 

01 

-6.4271E 

01 4.95426 

00 

— 5. 5438 E 

01 

ROW 4 











- 3* 0039E 

01 

-2.06266 

01 

1.7884E 

00 

1.0660E 

02 3 • 1 1 Ob 6 

oi 

-4.8 080 6" 

01 " 

-6.4176E 

00 

1.3715E 

02 

-3. 1614E 

01 

8.98136 

01 -9.57856 

00 

-2.6748E 

01 

ROW 5 











-4.19426 

00 

-6. 07 0 IE 

00 

5. 2 39 76 

01 

3. 11066 

01 3.5404E 

03 

— 4.6933 E 

02 

-1.98386 

02 

—9. 078 86 

02 

-1.86326 

02 

-6.03106 

02 -1.11906 

02 

-1.04156 

02 ' 

ROW 6 











2.09016 

01 

~TT0437E 

01 

— 6.8867E 

0 1 

-4.80806 

01 -4 .69336 

02 

7.9956E 

02 

9.72086 

01 

7.3279E 

02 

2.06826 

02 

6.89656 

02 1 . 8020E 

02 

2. 0033 E 

02 

ROW 7 











3.41216 

01 

1.1035E 

01 

-4.58026 

01 

-6.41766 

00 -1.98386 

02 

9.72086 

01 

3.47586 

02 

0.00006- 

01 

O.OOOOE- 

-01 

Q.OOOOE-Ol O.OOOOE- 

■01 

O.OOOOE- 

01 

ROW 8 











— 2. 4087E 

02 

-3. 9796E 

01 

— 1. 8365E 

01 

1.3715F 

02 -9.07886 

02 

7. 32796" 

02 

O.OOOOE- 

01 

4.06876 

03 

O.OOOOE 

-01 

O.OOOOE- 

-01 0.00006- 

■01 

O.OOOOE- 

01 

ROW 9 











4.33326 

01 

2.95386 

01 

— 3 • 2 17 56 

01 

-3. 1614E 

01 -1.86326 

02 

2.0682E 

02 

"o.boooe- 

01 

O.OOOOE- 

01 

2.82066 

02 

O.OOOOE- 

-01 O.OOOOE- 

01 

O.OOOOE- 

■or 

ROW 10 











-L. 92816 

02 

-1.0223E 

02 

-6.4271E 

01 

3.98136 

01 -6.0310E 

02 

6. 8965 E 

02 

O.OOOOE- 

01 

O.OOOOE- 

01 

O.OOOOE 

-01 

1.6917E 

03 O.OOOOE- 

•01 

0.00006- 

•01 

R0W1 1 











4.13466 

01 

3.61916 

01 

4.954 3E 

00 

-9.57866 

00 -1.11906 

02 

1.8020E 

02 

O.OOOOE- 

01 

0.00006- 

01 

O.OOOOE 

-01 

0.00006- 

-01 1.5570E 

02 

O.OOOOE- 

■01 “ 

ROW 12 











‘ -6.2867E 

01 

-6.09416 

01 

-5.5438E 

01 

-2.67486 

01 -1.04156" 

02 

2.00336 

02 

O.OOOOE- 

01 

O.OOOOE- 

01 

O.OOOOE 

-01 

O.OOOOE 

-01 O.OOOOE- 

•Oi 

2. 5553 E 

02 

COMPILE T IME= 

0.51 

SEC, EXECUTION 

T I ME = 

5.76 SEC, OBJECT CODE= 

416 
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MACH 0.9 


MATRIX ’D* 12 x 12 

THIS MATRIX IS NULL, 




MACH 0.9 


MATRIX ’K' 12 x 12 


ROW 

l 

COLUMN 

i 

1 • 739400D 

06 

ROW 

2 

COLUMN 

2 

2 .3348000 

06 

ROW 

3 

COLUMN 

3 

1.005800D 

07 

ROW 

4 

COLUMN 

4 

5.611200D 

06 

ROW 

5 

COLUMN 

5 

2.832599D 

08 

ROW 

6 

COLUMN 

6 

2 . 7 10200D 

07 

ROWS 

7 

THROUGH 12 

ARE ZERO. 




D3-WH4 


Y> 



MACH 0.9 


MA2RIX ‘Cjl* 12 X 12 
ROW 1 

9*546 If 03 5. 5596E 03 7.0879E 03 

2*78846 04 -3.6098E 04 7.0287E 03 

ROW 2 

•“ 4*47236 03 4.0766E 03 3.9739E 03 
7*32116 03 1.8704E 04 -3.2236E 03 


ROW 3 

' - 1.70746 03 2.8625E 03 5.45106 0 3 

4.0344E 03 -8.2189E 04 9.73076 03 

ROW 4 

~ 7.681 OE 02 6.7414E 02 1.20606 03 

-5.61766 03 -2.7928E 03 -6.15286 03 

ROW 5 

-7* 5 19 IE 03 2.8620E 03 1.6969E 04 
9.04236 04 -7.0943E 05 9*95196 04 

ROW 6 

-2.64596 03 -5.8633E 03 1.7116E 04 

-3.68826 04 2.3279E 05 -3.5577E 04 

ROW 7 

— 1.7955E 03 -7.86786 01 -7.07236 03 
-2.27526 04 1.5858E 05 -6.80656 03 

ROW 8 

2.5366E 03 1.2065E 03 9.1141E 03 

-i. 95016 04 2.672 IE 05 -4.49766 04 

ROW 9 

“ 2.1330E 03 -8.1875E 02 -5.2956E 02 
-1.9643E 04 7. 815 IE 04 -2.5 8556 03 

ROW 10 

-5.6653E 03 -6.5935E 03 4.8514E 03 

-1.6154E 04 4.6378E 04 -2.1548E 04 

ROW 11 

8.5478E 02 4.9323E 01 8.0694E 02 

-1.50506 04 5.9292E 04 -7.70076 03 


ROW 12 

-4.9659E 03 -5.5028E 03 -2.28445 03 
-5.06686 03 9.1340E 03 -4.5077E 03 


3.8111E 02 7.3005E 03 -4.1397E 04 

1.4390E 03 -4.5136E 03 1.7213E 04 


2.1366E 03 -7.59086 03 1.7285E 04 

2.5318E 04 -3.6786E 02 1.1730E 04 


3.02926 03 -1.523 IE 04 4.0480E~03 

7.66266 04 8.6543E 03 -5.7523E 04 


3. 0355 E 03 -1.0989E 04 2.05406 04 

8.0503 E 02 -6.7023E 02 -3.35266 03 


2.12806 04 -2.85456 03 -2.90016 05 
5.70826 05 2.3426E 03 -3.3593E 05 


8.18286 03 8.6830E 04 3.5704E 04 
1.4562E 05 9.5687E 03 1.01646 05 


4.5973E 03 -7.54946 02 9.0594E 04 

1.28756 05 -4.45546 03 7.3911E 04 


-2.0548E 03 1.3974E 04 5.07386 04 
1.46026 05 -1.6760E 04 1.3157E 05 


-4.001 2E 03 1.6280E 04 4.86626 04 

7.8550E 04 -4.38716 03 5.7091E 04 


4.0797E 03 4.00656 04 -1.29306 04 
5.0687E 04 2.4103E 03 1.84286 04 


-2.6278E 03 1.30126 04 2.72336 04 
5.67226 04 1.3131E 02 3.1074E 04 


-2.0469E 03 1.5074E 04 -4.18076 03 

-2.07006 04 2.72516 03 -8.39536 03 


'NO. ' P3-gBS5- 
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MACH 0.9 


MATRIX ’Cg 1 12 X 12 

ROW 1 

1.6705E 03 7.9264E 02 

-1.8203E 03 -6*575 2E 02 

ROW 2 

8. 88 ABE 02 7.4686E 02 

-3.5372E 02 -1.7343E 03 

ROW 3 

3.4490E 02 4*43348 0? 

-1.4008E 01 5.5340E 03 

ROW 4 

-1.6831E 02 9* 2 80 IE Oi 

4.8746E 02 2.0345E 03 

ROW 5 

9.9021E 02 -4.5048E 02 
-5.6282E 03 3.6540E 04 

ROW 6 

—6* 521 IE 02 3.1337E 02 

2.1516E 03 -8.3494E 03 

R OW 7 

1.6465E 02 8* 561 58 01 

2.5817E 03 - 1* 0264E 04 

ROW 8 

-2.0028E 03 -4*751 OE 02 
1.5568E C3 2. 1085E 04 

ROW 9 

1.9987E 02 2.7050E 02 

1*53176 03 -5.3167E 03 

ROW 10 

-1.2489E 03 -5 *6559E 02 
1.1489E 03 -1.5265E 03 

ROW 11 

3*51 47E 02 3.8958E 02 

l* 0868E 03 -3* 1943E 03 

ROW 12 

- 2 • 8 02 9 E 02 -1.1244E 02 
3.3752E 02 -1.6897E 03 


5* 41 3 0E 01 -1.4970E Oi 
-7*9383E 02 -2*2632E 03 


3*90298 02 3.8410E 02 

-8.0198E 01 -2* 3859E 03 


7* 991 OE 02 2 *55198 02 

— 5*9 770E 02 4*7989E 03 


2.7861E 02 5 • 0808 E 02 

7.2578E 02 1.7216E 03 


3 * 061 2E 02 -2*04438 03 
-5.4079E 03 2.9621E 04 


-1.1184E 03 1.3945E 02 

1.2779E 03 -2 • 7962E 03 


7.2780E 01 2.4389E 02 

3.7307E 02 -3.3874E 03 


-5.1813E 02 1.1778E 03 

2 *79 16E 03 -8.8453E 03 


-1.8201E 02 1.1540E 01 

1*0223= 03 -5.1715E 03 


-9*49956 02 -1.3743E 01 
1 • 2939E 03 1.6583E 04 


5.1264E 01 7 *96498 01 

7.3777E 02 -3.1029E 03 


- 4* 963 IE 02 -2.4835E 02 
2.6709E 02 1*5070E 03 


1.3304E 03 -2.7772E 02 
-2*9289E 02 -2.5975E 03 


9.1696E 02 —2 *9966 E 03 
-6 *2382E 02 -1.7279E 03 


2.3679E 03-1*65406 03*" 
-9.9402E 02 2. 5925E 03 


l • 1302E 02 -1.9468E 03 
-1.6800E 02 -3.3669E 02 


1.1791E 04 1.2904E 04 

3.1472E 02 1.9463E 04 


-1.3139E 04 1.1394E 04 

-7 *87 00 E 02 -3* 1 175E 03 


-i*5685E 03 -5* 3703 E 03 
2 *7 486E 02 -4.8496E 03 


-6*51 23E 03 1.0508E 01 

8 *95058 02 -9.5699E 03 


-2 *4154E 03 -l* 7568E 03 
1.3303E 02 -3* 7402E 03 


—7*40488 03 6.9657E 03 

2 *0522E 02 -4*2349E 03 


— 1 *9351E 03 -B.4662E 01 
4.5643E 02 -i*5648E 03 


-2 • 9284E 03 3.9663E 03 

-2 *99316 01 6.5071E 03 



NO. D3-888U- 
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mach 0.9 


MATRIX 'C 3 * 12 x l2 

ROM 1 

7*60106 00-3*57116 00 4.3557c Oi 

- 5*65676-01 -4.36976 02 6*56066 00 

ROM 2 

3.48166 00-2.75366 00 2.3084E 01 

-4.33226 00 -1.42056 02 -5.1231E 00 

ROM 3 

-6.1740E 00-3.02116-01 -6.4186E 00 
1.64816 01 1.9483E-01 2.95196 01 


ROM 4 

-2.18376 00 -3.4294E 00 -2*08956 00 
- 7.88026 00 1.3042E 02 -2.1445E 00 

ROM 5 

-2.51776 01 -7. 2157E-01 -5.8300E 01 
1.29636 02 -6.5390E 02 1.8493E 02 


ROM 6 

1.79396 01 1.8220E 01 5.0111E 01 

-1.9063E 01 5.1092E 02 -1.3536E 02 

ROW 7 

3.5358E 00 -3.8162E 00 1.6843E 01 

-4.85636 02 -3.0861E 01 -3.0157E 01 

ROW 8 

-7.8735E— 01 1.5306E 00 -1.15156 01 

-2.37686 00 9.241 8E 02 -6.2697E 00 

ROW 9 

”8.68896 00 -4.57016 00 2.4964E 01 

-3.62516 01 -4.2471E 01-4.15026 02 

ROW 10 

2.3839E-01 2*69906 00 -6.7166E 00 

-9.72306-01 2.73846 02 -3.26166 00 


ROW 11 

8.33986 00 1.57156 00 2.8863E 01 

-4.2305E OC -5.33896 01 -3.4549E 01 

ROW 12 

7.6270E-01 2.40186 00 -2.4082E 00 

-1.91 60E-01 3.65606 01 -7.1538E-01 


~ ■■ ■ 






1.12876 

01 

1 .75 356 

02 

-2.7407E 

02 

-5.15036 

02 

2 .65476 

00 

-4.60986 

02 

9.33266 

00 

1.07476 

02 

-1.51586 

02 

-2.63876 

02 

2.0084E 

00 

-3.27166 

02 

6.9502E 

00 

-5.12 826 

01 

-8. 0458 E 

00” 

-1.07366 

02 

3.35366 

01 

-2.60966 

02 

4.02306-01 

-2.0547E 

01 

3.2390E 

01 

1.05466 

02 

1.8164E 

01 

-1.16266 

02 

-7.9895E 

00 

- 4 .33296 

02 

2.22386 

02 

-7.12816 

02 

i .7 4026 

02 

-3.93916 

02 

2.58036 

01 

3.1127E 

02 

-4 .1082E 

02 

7.66456 

02 

- 2 • 3061 E 

02 

7.0130E 

02 

-6.43726 

00 

1.1436E 

02 

-2.T478E 

01 " 

-2.5608E 

01 

-4.58166 

oo 

-1.53796 

01 

-7.61266* 

-01 

-1.60446 

01 

6.42436 

01“ 

4.68136 

02 

-1.0459E 

01 

1.78286 

02 

-2.23766 

00 

l .48916 

02 

-1.0665E 

02 “ 

-3.50506 

01 

-3 *07316 

01 

-2.13116 

01 

-1.11336 

00 

-1 .4352 E 

01 

3 .79536 

01 

7.0704E 

02 

-8.84036 

00 

3.00166 

02 

1.1887E 

01 

1.47746 

02 

-1.84686 

02 

-4 • 8901 E 

Oi 

-3 .73776 

02 

-2.97986 

01 

-3. 17656 

-OE 

-4.1338E 

00 

1.00876 

01 

9.848 IE 

OIL 

-2.69886 

00 

3.24556 

02 





w .. . 
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MACH 0.9 



MATRIX ’D^ 


12 X 12 


ROW l 

-3.8963E-01 -2.6581E-02 
l* 0733E 00 — 1 • 8130E 01 

ROW 2 

2.8259E-02 7.3967E-02 

6.3829E-01 -5* 7048E 00 

ROW 3 

" -9.6866E-02 -3.9988E-03 
-I.7867E 00 1.0673E 01 

ROW 4 

8*469 5E-02 1.7019E-02 

7.1957E-02 2.5413E-01 

ROW 5 

-2* 345 IE 00 -6 • 4414E-0 1 
-8.9257E 00 7.3517E 00 

ROW 6 

-2.978 4E-01 2.2923E-01 

3.7764E 00 -1.8142E 01 

ROW 7 

7.9961E-01 I. 231 OE-Oi 
6* 841 4E-01 9.6328E 00 

ROW 8 

* 1.2416E-C1 3.0556E-01 

1 *6544E 00 -6* 065 OE 00 

ROW 9 

3.5254E-01 -2.4632E-02 
1.4881E 00 -4.9205E 00 

ROW 10 

-2.8069E-02 1.3026E-02 

9.4530E-01 -2.2390E 00 

ROW 11 

1.7776E-01 1*015 IE- 01 

l . 27 CAE 00 -l • 78 1 3E 00 

ROW 1 2 

-1.2377E-02 -1.4488E-02 
2.3903E-0L -4.2867E-01 


7* 4763E-02 5.8770E-G1 

2.4185E 00 -1.7237E 01 


3* 8757E-01 2.1889E-01 

7.0565E-01 -3*4562E 00 


-1.3268E 00 —2 *6053 E— 01 
-4.5425E-01 6.9562E 00 


2.0733E-01 -4.9515E-02 
2.9475E-01 8 • 1360E— 0 1 


-8.4280E 00 -2 *3052E— 01 
3* 4887E 00 -2. 1637E 01 


1.7037E 00 7.2501E-01 

1.4225E 00 ~1*0432E 01 


1.7484E 00 -3.3568E-01 
2.2884E 00 1.5128E 01 


1 *452 IE 00 5.6885E-01 

5*278 6E -01 i • 1 074E 00 


1 • 5462E 00 - 1 • 0479 E-OI 
-9 • 1957E-0 1 -5.6152E-01 


6* 5387E-0 1 1.4849E-01 

6* l 045E-03 9.8376E-01 


9.8143E-01 9.2114E-02 

-2*10 70E-01 1 • 1444 E 00 


1*340 3E -01 -S.9169E-05 
-5.8060E-02 3.9188E-01 


-1.7933E 

00 

-5.3592E 

00 

1.3049E 

00 

-1.0905E 

01 

6.6017E-01 

-3.4866E 

00 

-2 *4847E-01 

-4.3527E- 

-01 

-3.0646E 

00 

l *0054E 01 

9* 1848E* 

-01 

2.I138E 

00 

5.0360E- 

-01 

-5.8009E- 

-01 

— 3 • l 209E— 01 

8* 0377E— 01 

-2.6573E 

01 

4.7544E 

01 

1.1868E 

01 

-3.7099E 

01 

3.9927E 

00 

-1.6180E 

01 

-1 • 8976E 

00 

“3.8878E 

00 

6.8646E 

00 

-6*2 642 E 

00 ' 

-2.4364E 

00 

1.5876E 

01 

5.4383E 

00 

- i • 0869 E 

01 

-2.1719E 

00 

6.6897E 

00 

4.0278E 

00 

-7.7825E 

00 ■- 

-1.8100E 

00 

3.0364E 

00 

2.4038E 

00 

-4.4033E 

00 

-I *1272E 

00 

2.0941E 

00 

2 *7 845E 

00 

-6* 5899E 

00 

-I *1 2&7E 

00 

2.9155E 

00 

4.4168E 

“01 

-7.2738E-01 

-2.7893E 

-01 

4.4718E 

-01 


53-8884 
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MACH 0.9 

MATRIX *ch * . ~ ~ 

l X 12 

ROW 1 

3* 00006-03 3.0000E-03 3.0000E-03 3.0000E-03 3.0000E-03 3.0000E-03 

3.0000E-03 3. OO0OE-O3 3.0000E-03 3.0000E-03 3.0000E-03 3.0000E-03 




MACH 0.9 


matrix 'd 2 

f 

12 X 12 










ROW 1 

1.4513E 

01 

3.2915E 

00 

1.5738E 

01 

-1.4290E 

Cl 

9.9497E 

01 

2.6848E 

01 

-1.0380E 

01 

4.4485E 

02 

-6.8389E 

01 

4* 9389 E 

02 

-5.6169E 

01 

3.7616E 

02 

ROW 2 

4. 3993 E 

00 

4* 6936E * 

-01 

-1.3576E 

01 

- 1 * 041 IE 

01 

3.9756E-01 

1.4967E 

02 

-2.3812E 

01 

3.9244E 

02 

-3.7696E 

01 

3 • 4658 E 

02 

-1 .0769E 

01 

2.1586E 

02 

ROW 3 

-8.7769E 

00 

-2. 7055E 

00 

3.0759E 

01 

1.5446F 

01 

3.6979E 

01 

-3.2677E 

0 2 

5.94266 

01 

-5* 7547E 

02 

4.3259E 

01 

-5. 1664E 

02 

4.7826E 

00 

-3.3679E 

02 

ROW 4 

-1.5242E 

00 

-9.8026E- 

-01 

-6.3062E 

00 

3.0322E- 

-01 

-7.9325E 

00 

2. 9946 E 

01 " 

-5.2177E 

00 

2.9415E 

01 

4.9833E 

00 

8.8156E 

00 

3. 161 IE 

00 

-3.0872E- 

■01 

ROW 5 

-8.63 77E 

00 

-5. 0649E 

00 

3 • 17 1 IE 

02 

8.5604E 

01 

6.'+598E 

02 

-2.5106E 

03 

4.2669E 

02 

-3.5323E 

03 

1 .70248 

02 

-2.7149E 

03 

-1.5237E 

02 

-1.4452E 

03 

ROW 6 

2. 1308E 

01 

-3.634 2E 

00 

-4. Ill 4E 

01 

-3 » 1982E 

01 

-1.3653E 

02 

4.8447E 

02 

-1.1331E 

02 

5.6909E 

02 

-5.3253E 

01 

4.0959E 

02 

4.9099E 

01 

2.2237E 

02 

ROW 7 

-4.2653E 

00 

5* 2496E 

00 

-9 .22906 

0 1 

— l • 1134E 

01 

-2.0177E 

02 

6. 13076 

02 

-7.8175E 

01 

8.6L53E 

02 

-1.2240E 

00 

6.66696 

02 

1.5169E 

01 

3.4813E 

02 

ROW 8 

" 6.7823E 

00 

-7.6610E 

00 

-4.4903E 

01 

-2.8469E 

01 

-1 .5105E 

02 

4.1 634E 

02 

-6. 442 IE 

01 

4.4428E 

02 

-4.7042E 

01 

2.3695E 

02 

4.5309E 

01 

-2.1259E 

01 

ROW 9 

4.6182E 

00 

6.1454E 

00 

-5 .9905E 

01 

-1.0686E 

01 

-9 .3745 E 

01 

4.4013E 

02 

-7.5779C 

01 

7.8273E 

02 

-2.0423E 

01 

6.7461E 

02 

2 .24676 

01 

4.2890E 

02 

ROW 10 

6.6302E 

00 

l • 2 60 9E 

00 

-2.0685E 

01 

-9. 3781 C 

00 

-6.6337E 

01 

1.7413E 

02 

-3.7415E 

01 

2. 007 IE 

02 

-1.5159E 

01 

1.0172E 

02 

2.5506E 

01 

3. 13616 

01 

ROW 11 
4.3949E 

00 

2.0971E 

-01 

-3.6235E 

01 

-1.2319E 

01 

-7.6890E 

01 

3. 189 IE 

02 

-5.3938E 

01 

4 » 0033E 

02 

-2.2912E 

01 

3.3663E 

02 

9.4659E 

00 

2.0770E 

02 

ROW 12 

1.7463E 

00 

8. 840 46 

-01 

-4. 29756 

00 

-9.9553E 

-01 

-1.1254E 

01 

3.0617E 

01 

-9.6478E 

00 

4.6 58 OE 

01 

-1.7337E 

00 

1.9076C 

01 

6.7206E 

00 

8.5784E 

00 


JBF49JF * fA/£2 NO* D3-8884 
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mach 0.9 


MATRIX , d 2 * 1 X 12 

ROM 1 

9.0000E-03 9* 0000E-03 9.0000E-03 9.0000E-03 9.0000E-03 9.0000E-03 

9.0000E-03 9.0000E-03 9.0000E-03 9.0000E-03 9.0000E-03 9.0000E-03 


NO. D 3-8884 




MATRIX 'D 3 ' 12 X 12 

ROW 1 

-4.1001E 01 -8.5356E 00 -t.3365E 02 
-3.7A58E 01 -8.8A6AE 02 2.3A31E 02 

ROW 2 

-1.7746E 01 4.1547E 00 -2.9263E 00 

6* 7257E 01 -1.1862E 03 1.5974E 02 

ROW 3 

4.X318E 01 5.2318E 00 -4*3543E 01 

-2.0159E 02 1« 472 BE 03 -1.7634E 02 

ROW 4 

5.6176E 00 7.1805E 00 2.4934E 01 

2.4S60E 01 -8*095 IE 01 -2.0626E 01 

ROW 5 

1.4432E 02 -2.5445E 01 -7.8839E 02 
-1.7170E 03 l * 1600E 04 -9.6005E 02 

ROW 6 

— 1.1306E 02 1*620 IE 01 9. 8463E 01 

5.3704E 02 -1.3334E 03 3.0068E 02 

ROW 7 

-1.4146E 01 -8.1473E 00 2.4280E 02 

3.4671E 02 -3.0370E 03 1.941 8E 01 

ROW 8 

-6* 231 4E 01 2.6527E 01 1.6262E 02 

2.8194E 02 -2.0198E 03 2.7591E 02 

ROW 9 

-3.3757E 01 -I.O047E 01 1.3141E 02 

3.1736E 02 -2.2366E 03 1.2327E 02 

ROW 10 

-5.0625E 01 -1.3906E 01 7.3491E 01 

1.7606E 02 -1.0324E 03 1.0770E 02 

ROW 11 

-2.8068E 01 8.2371E 00 7.4442E 01 

2.1927E 02 -8.9913E 02 1.2071E 02 

ROW 1 2 

-1.2839E 01 -5.8845E 00 1.4336E 01 

4* 71 96E 01 —2 • 5649E 02 1.7967E 01 


MACH 0.9 


3.0157E Oi -5.7912E 02 4*9503E 02 

-1.2067E 03 2.4141E 02 -1.0393E 03 


3.6125E 01 -1.8744E 02 -2.5595E 02 
-1.09UE 03 7.0214E 01 -7.2060E 02 


-5.9188E 01 -6*2 7 69E 00 8. 3 985 E 02 

1.3894E 03 -3.8601E 01 9.3225E 02 


1.6915E 00 3.5982F 01 -1.2B34E 02 

7.3255E 00 -4.4439E 01 3.2904E 01 


-4.1060E 02 —1.3 081 E 03 8.0301E 03 
9.0615E 03 3.8217E 02 4.8568E 03 


1.6680E 02 5.8640E 02 -1.7514E 03 

-7.7404E 02 -2.2462E 02 -1.6100E 02 


5.4049E 01 3.7930E 02 -1.9599E 03 

-2*4201 E 03 3 • 1486E 01 -1.3323E 03 


1.5054E 02 6.0249E 02 -1.7249E 03 

-1.2339E 03 -1.4132E 02 -1.8575E 00 


4.9288E 01 9.7467E 01 -1.2630E 03 

-2r0011E 03 -9.6342E 01 -1.3447E 03 


5.1182E 01 2.3122E 02 -7.0332E 02 

-6.2517E 02 — 8 *3946E 01 -3.0074E 02 


5.8097E 01 1.7330E 02 -9.7675E 02 

-7.7190E 02 -6 • 05 88 E 00 -4.5589E 02 


6. 1077E 00 2 *7 504E 01 -1.1222F 02 

-1.2902E 02 -2.3873E 01 -7.7980E 01 



• 

o 
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MACH 0.9 

. matrix m 3 ‘ t x l2 _ 

ROW 1 

1.5000E-02 1.5000E-02 1.5000E-02 1.5000E-02 1.5000E-02 1.5000E-02 

!• 5000E-02 1 •50Q0E-02 1.5000E-02 1.5000E-02 1.5000E-02 1.5000E-02 


NO. p 3 - 8884 
[sect PAGE 3 ° 





MATR 1 X 

» 

12 X 12 


ROW 1 
2* 4698E 

01 

6.2633E 

00 

8*49858 

01 

6.3286E 

02 

ROM 2 

1.05476 

01 

-8.5X82E 

00 

-4.1052E 

01 

9.8037E 

02 

ROW 3 

-4.50646 

0 1 

-8.3073E 

00 

1 * 8080E 

02 

-1.2040E 

03 

ROW 4 

-5.6863E 

00 

-1.0229E 

01 

-3.0626E 

01 

1.4945E 

01 

ROW 5 

-1.9291E 

0 2 

2.066 7E 

01 

1.6618E 

03 

-1*041 5E 

04 

ROW 6 

1.5449E 

C2 

9.2314E 

00 

-6.0236E 

02 

1.6778E 

03 

ROW 7 

2.3423E 

01 

5.4258E 

00 

-3.7359E 

02 

2.7116E 

03 

ROW 8 

9.6640E 

01 

-1.2 599E 

01 

— 2* 8S88E 

02 

2.3158E 

03 

ROW 9 

3.8510E 

01 

7 *63096 

00 

-3.2352E 

02 

1.7796E 

03 

ROW 10 

7.3 539E 

01 

2.6220E 

01 

-1.9061E 

02 

1.2503E 

03 

ROW 11 

3.5752E 

01 

-5.7490E 

00 

-2.1832E 

02 

6.9490E 

02 

ROW 12 

1.8074E 

01 

8.5987E 

00 

-5.2434E 

01 

3.1915E 

02 


1.52496 02 -9.8617E 00 
- 1 . 993 66 02 9.7672E 02 


2.3376E 01 -2.8108E 01 
-1.5019E 02 9.1038E 0? 


2.4608E 01 5*6734E Oi 

U7815E 02 -l • 1734 E 03 


- 2.735 OE 01 -5* 8881 E 00 
1.7770E 01 -7.9392E 01 


6.2636E 02 4.2987E 02 

1.0790E 03 -8.2129E 03 


-9*5 128E 01 -1.9228E 02 
-4.2269E 02 1.1134E 03 


-1.9617E 02 -5.3 145 E 01 
1.0205E 01 2.1837E 03 


-1.5650E 02 -1.7139E 02 
-3.4044E 02 1.6162E 03 


-1.0128E 02 -5.3275E 01 
-1*631 2E 02 1.6300E 03 


-6.3283E 01 -5.5990E 01 
-1.4412E 02 8.4979E 02 


-5.7555E 01 —6* 1898E 01 
-1.4062E 02 6 • 1838 E 02 


-1.0649E 01 -6* 2147E 00 
-2*68236 01 1.7576E 02 




MACH 0.9 


6 *68 95E 

02 * 

-7* 2666E 

02 

-2 .44046 

02 

9 * 222 IE 

02 

2 • 6 0 3 2 E 

02 

6.9608E 

01 

-7.4982E 

01 

6.2376E 

02 

2.2151E 

00 

-6.B127E- 

02 

4*95776 

01 

-7.90616 

02 

-6 .0655E 

01 

1 *51946 

02 

5*77856 

01 

-9.92476 

01 

l .03186 

03 

-6.9596E 

03 

-2 .8032E 

02 

-4.3607E 

03 

-5 .86 1 3 E 

02 

1.7564E 

03 ~ 

1 .6604E 

02 

3*12616 

02 

-2.7970E 

02 

1*66626 

03 

-5.3675E 

01 

1.2245E 

03 

-5 *76966 

02 

1.71716 

03 

9.0735E 

01 

2.2684E 

02 

-7.0402E 

01 

1.0660E 

03 " 

1 .1744E 

02 

1*1267E 

03 

-1.5952E 

02 

6 . 37 36 E 

02 

5.9075E 

01 

5.2484E 

02 

-1 .6771E 

02 

8.75246 

02 

-2.0140E 

01 

3.4229E 

02 

1 1.3903E 

00 

8.0755E 

01 

1 1.9267E 

01 

1 • 2423E 

02 


| ; ]Q, D3-8884 
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MACH 0,9 

MATRIX f d4' 1 X 12 

ROW 1 

2* I'OOOE— 02 2*1 000E-02 2.1000E-02 2.1000E-02 2.1000E-02 2.1000E-02 

2.1000E-02 2.1000E-02 2.1000E-02 2.1000E-02 2.1000E-02 2.1000E-02 




MACII 0.6 


MATRIX •CjL* 8X8 


ROW 1 

2.O201E 04 1 * 05C1E O4-8.8509F 02 -4.0262E 03 -1.4376E 04 2.1562E 04 

5. 11546 03 6* 9209E 03 


ROW 2 

1.04R9E 04 8.3595E 03 2.5310F 03 -8. 136*5? 02 -1 . 2*556?- 04 1.9407E 04 

2.6551E 03 7. 9292? 03 


R DW 3 

1.R045E C3 5.0609F 03 7, 8430E 03 2.7232C 0? 5./608E 03 4.4630E 03 

-6.8414E 02 - 1 • 2003? 03 

R -l.5160E 03 6.5938F 02 4. 1986E 03 6.2190E 03 -4.7119E 00 -9.6726F 03 

-2.7646E 03 5.5280E 02 


^°6.4537E 03 1.9413E 03 9.0176E 03 -5.1823F 03 7.9361E 04 -6.69696 04 

-4.7343F 03 -4.4895E 04 __ 

.^^ 7 ^5 e ~04^5i-7836E~03 — 6v0567E- 03 - 6* 4927E-02 —6' v8556E- ** ls-i998 6 05" 

1.6866E C4 6.8343E 03 

- R °1.546 7 E 02- -1.5942F 02 -1.0369E 03 -1.4290E 03 -5.7554E 03 1.09336 04- 

9. 14 59 E C3 4 . 448 IE 03 


-l.? n 04E C4 -6.791 8 E 0? -2.3800E 03 4. 6653E C3 -1.9241E 04 1.02346 04 

-3.5203E 02 5.008 1£ 04 


I NO. D3-8884 
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MACH 0.6 


MATRIX ’Cg’ 8X 8 _ 

ROW 1 

- 1.5608E- CS— 6;4fc6ir 0? 9.5133E 02 3. 1098 E 02 5.79Mf 03 -7. 2531 F 03- 

-9.4633E C2 -3.9740E 03 _ _ 

ROW 2 

7i36 31E“02" ~5r6'i ?6E ~0^ 7w 51 CAE 02 — 4v 0109E 02 ~2 . 9fr54f-03 “4v 8371 O ' 03 

-5. 30 CAE 02 -2* 51C7E 03 . 

ROW 3 

-1.7231E 02 2m2t22F 02 ; 8.9C28F 02 “ 4. 2083E 02 iv7420E 03 -2.7360E 03 

2.1228E 02 - 1* 2695E 03 

ROW 4 

-2-2652EC2 -3 *7073 ti'Ol *-5.9045E Ot 3.3301P 02 — 1 *22156 C3 3 *95746 02 
5* 563 IE 02 1-1 8146 03 _ 

ROW 5 

- j . 02 0 1 £-03 -— IvCAe ie-03 — 1. 59736 03 — l . 60 136~02 -0.26166 03 -5. 14886- 03- 

1.1920E 03 -7.74286 03 __ 

ROW 6 

l*S271E 03 - 1. A735E 03 4.0551E 0? -6.6469E 02 -4.3669F 03 ■ 9* 3129 E 03 

-2.287AE 03 1.31926 OA 

R OW 7 

’7.24-58F 02 4.1441E 02 -4.0210C O? -4.2103E 0? -2.0537E C3 1.5752E 03 

-6. A3 99 E 02 -R. 92566 01 

ROW 8 

-3.4428F C2 -1.3M4E 02 -8.P177E 0? — 1 . 52E5F 0? -3 . 7766E- 03 5. 0(106 B O* 

9. 849 IE 01 1 • 4307*F 04 



NO. 

D3-8884 
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8X8 



MACH 0.6 


M ATP IX *C^ ? 

ROW 1 

6.76426 00 -4.13946-01 
3.24976 CC -3.76386 C2 



3.4064F" 01 9. 66796 00 1 .39936- 02 — 2.16 26F- 02 


ROW 2 

— 3. 38 896~O0 — A . 95 1996-01- 1.91486-01 9. C690t~00 8* 735 2 6-01 l-r29&<re-02-- 

-3.74?6E 00 -2.047SE 02 


ROW 3 

-5.46546 00 1. 1496>6- 00 -7.97886 00 6*01046 00 -5.24296 01 -3.33006 00 

2.5414E 01 -5.45856 01 


ROW 4 

-2.0698E 00 -2.656356 OO -2.4685F OO 5.3893F-01 -2.21436 01 3.11696 01 

-3.24126-01 8.96136 01 

ROW 5 

2. 39 20E 01 2.41716 00 -6.55166 01 - 1.31 706 01 - -4 .1 5 0 1 6 - 02 — - 2 . 8 404 E~0? - 

1. 62 1 IE 02 -3.59C7F 02 


ROW 6 

1.5433F 01 1.08246 01 6.5762E 01 3.24586 Cl 3.14156 02 -4.41566 C7 

- 1 • 3^ 92E O? 4.33826 02 

ROW 7 

- 9.-94C86 00 -3.75976 00 2 . 87886 01 -1 . 4385F OO “ 1 .6 807E 02 -1.33906 02 
-4.4658E 02 -5.9373E 01 

RCW 8 

-4. 88 59F-C1- -4.94576-01 -3.912 5E 00 -5. 77246-01 -8.76666 00 - 2.60316 01 

— 3.8031E 00 5.2264E 0? 


MATRIX , D 1 * 


ROW 1 

3 v8 8 14F-02 " 4.4 4266-0? -8.2333 6-0? — 9* 1RC6F-03 —3.54506—01 1. 0764 E— 01 - 

2.3468 E-C2 2.37956-01 

ROW 2 

3. 86856-02 2.63446-0? 1.8294E-0? -1. 94796-02 5.C-892E-02 1.1487F-01 

-3.1710F-C2 3.581716-01 


ROW 3 

■ 4.5372E-02 - 3 . 499 6.9160F-0? -1.71166-0? 1.86066-01 2.5142E-01 

-2.0123F-C2 1 • 7698 F— 01 

ROW 4 

-2.75466-0? -8. 8435F-03 8. 149 3F-03 - 1 . 08C 1 6-0? -2*37666-03 -8.55056-02 

1 • ?0 4 1 F— 0? -1.4 34 76-01 


ROW 5 

4 . 5 3 4 5 F. - 0 2 -4.265?£-''? -3. 17036-01 1. 1621F-C1 -5.5635E-C1 -4.7284F-^1 

1.97856-01 -1.346CF 00 

ROW 6 

1.27 66 F-^l 1.45426-0? -1.28676-01 -1.5990F-01 -3.19966-01 2.1220F 00 

-4.8409F-02 2.30876 CO 
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MACH 0.6 


ROW 7 

2.0771E-02 . 1.3274E-02 3.5358F-02 -5. 36 53E-C2 3 . 3928 E-02 

-4.4746F-C2 7.3141 E— C 1 


ROW 8 ' ^ 

-l.<T?50f-C2 -l.~2903F- 0?- 9.5901F-03 — ii2150E-02 — 2v€163E-03 

- 3 • 2 794E-03 4.7058F-0? 


MATRIX , d 1 * 


1X8 


ROW l • ■ ‘ M 

- 3.0C00E-C3 3.0000^-03 ~ 2 .0009F-03 3. OOCOF-03 ‘3 .00008-03 

3. 0O09E-C3 3.000CE-C3 


MATRIX *02’ 


8X8 


ROW 1 A _ 

1 *4591 F 00 -7.703CE-^>1 -2.3938E 00 -1.1741F 00 -6.6840F- 00 

6.1977E-C2 4.15575: 01 


-1.40I2E 00 -5.122 IE— 0 1 -4.2=5 14E 0^ 3. 1087E-01 -1.4409E 01 

1.9281F 00 1.5056E Cl 


ROW 3 

-2.0573E 00 -3.51971F-01 1.9999E On 1 . 1976F 00 8.8318E-01 - 

6.0818F-01 -3.1990E Cl 


ROW 4 

l.O^AOF 00 2.6345C-01 3.7<>77F-01 -3.5128F-01 2.6110F 00 

—6. 5304F-01 4.0047'F-ni 


ROW 5 

-3.5245F DO -1.1603# 00 3.05R8F. Cl -1.6684F OO 8.4048E 01 

-1.0672E Cl -1.0 539'IF 02 


ROW 6 

-5.2233E CO 7.641 61E-01 -6.91 20F 00 5.4C90E 00 -*2.9C6^E 01 

5.0433E 00 -6.7264E 00 

ROW 7 
-A . 22036— or 
3.16C5F 00 

ROW P 

7*8461 E— 0 1 
1. ?943*=- n l 

MATRIX , d2 f 1X8 


- 5i 1 21 2F-01 -3. F685F 00 1 i 1596F 00 -2*62506 01 

3 .980 5 E 01 


4.1376E-01 -1.6686F-02 5*4282E-P1 1.4314E 00 

-5.8504IF CO 


ROW 1 ■ 

O.^^OOF-O? 9.00006-^3 9.0000E-03 9.0OC0F-O3 9.0000E-03 

q.nnonp-03 9.noOOE- r ‘3 
MATRIX 8X8 

ROW 1 

9.1101F 00 7.464PE 00 -1.1236E 01 6.7142F 0^ -3.6324F 01 

1 H3-B884 Page 


• 6 261F— 01 


v6662 E—OI 


3.0000F— 0> 


3 .3461 F 01 

1.2 857 E 01" 


2.3747F 01 


1. 1867 F— 01 - 


7.5624E 01 


3. 1609 E 01 


1.7220E 01 


8.6161F 00 


9. OOOOF-93 


- 1 • 0097 E 02 

4? 



MACH 0.6 — • 


5.4919E CO -6*75026 Cl 


pdw 2 

7.0R87F 00 4.62266 00 1. 2940F 01 “l.2379fr 00 4.56806 01 -3.5677E 01- 

-8.1043F CO -1.6385E 01 


pnu 3 

1.1201E 01 9.88396-02 5.701 06 00 -6.73366 CO 5.0195E 01 7.69446 01- 

-8.764 7F 00 9. 19846 01 


ROW 4 

-6il335E 0& -Ii82«8E 00 - 8. 88668-01 t. 86538 00 -1.11088 01-8.36888-00— 
3, 7001 E CO -1.39818 01 

ROW 1 ? 

2.6367E 01 -4.91336 CO -1.0021E 02 6.5014E CO -2.21816 02 1.9556E 02 * 

3.7737E Cl 1.81676 C2 ■ 


R 2.89296 01“ 1.9345F 00 2.2368E 01 -3 . C390E Cl 1 . 1215E 02 2.50526 02 
-2.91C7E 01 2.0766E 02 


ROW 7 

— 9; 226 IE-01 - 1.3 774E 00 — 2.4874F 01~-6 rC340E- 00 - * 6 . 8356E Cl -1 . 55 1 7 F -01— 

-1.1*966 01 -2.3397F Cl 

ROW 8 

-8.53388 00 -3.30928 CO 2.88678 00 -3.C591E 00 -2.23908 00 3.97068 01 

-1.3037E ^0 2.0643F 01 ^ 

MATRIX 'd^* 1 X R 

p OW 1 

i. 5000E-02 1 • 50006-02 1.50096-0? 1.50C06-0? 1.5000E-C2 1.50COF-02 

1.5 r O r 'E-02 1 • 50CCE -02 

MATRIX ' d 4* 8 X8 


ROW 1 

-1.2524 F 01 - 6 . 9440 E OC 1 . 8970F 01 -8.7719E 00 5.6766F 01 1.1864F^2 

-7. 84986 Of 9.9696F 01 


ROW 2 

— 7. 1 ?8 f 'E 00 -5.5393F n 0 — 1.61736 01 6.18196— n l —6.351 IF Cl 5.4080E 01 


9.08066 

00 

3*36306 

Cl 








ROW 3 

-l. 4' 996 

Cl- 

-4 .70?o£- 

►01 

-1 . 4005E 

01- 8 .88426 

CO 

-6.66866 

01 

-8.51626 

01 

1.44996 

ot 

-1 . 1 190 E 

02 








ROW 4 

7.7479F 

00 

2.°785F 

0 ^ 

1. R585E 

0" -2.23456 

no 

1.6151F 

01 

3. 54O0F 

90 

-5.21616 

00 

1,66656 

01 






— 


ROW 5 

-4.08656 

01 

. 35466 

09 

1.14106 

0? -1. 1378E 

oo 

1 .92746 

02 

-2.93156 

02 

-2.82 03 E 

ni 

-3.41736 

C ? 




- 
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MACH 0.6 


ROW 6 

-3.0^770 01 -3*94476 00 -4.1630E 01 3*7872C 01 -1.6303E 02 -2.7053E 02 

4.0927E 01 - 2 • 1 B 5 76 0? 

ROW 7 

3,47 2?F”CO -i.3549>£ -9*03 546' 01 6*090 1£ OO - 6v7695E " Oi - 3 , 5904E 01 

1.13096 01 6.1649E 01 


ROW 8 

6.3647F 00 2.7?19»6 00 -3.6583E 09 3v811 6E 00 4.2R11E CO -4.7258E 01 

1.3847F 00 - 2 • 080 IE 01 


MATRIX 'd^ 1 1. X 8 


ROW 1 

?• 1003F-C2 2 • 10006— 0? ?i 1000F-C2 2. 10000-02 2.100CE-02 2,10006-02 

2.1000P-02 2. 1PC06-C2 
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MACH 1.2 


MATRIX *C 1 1 8x8 


2.6C763E 

03 

—3*633956 

02 

3.83336F 

03 

-4.29293E 

02 

9*039756 

03 

1.368C9E 

04 

3, 86283E 

04 

-2*517416 

04 







l. 666846 

C4 

1.45C29E 

03 

4.17266E 

01 

2*104346 

03 

1.10184E 

02 

7* 31426E 

02 

2*889338 

04 

-1*495766 

04 









7 • 378096 

02 

-4. 5668 CE 

01 

2*780776 

03 

1.11872E 

03 

-2.168756 

03 

6*864886 

03 

2*558108 

04 

-1*32 4826 

04 







1*567256 

03 

-2*02867 € 

02 

6.170166 

02 

_ l .391876 

03 

1 .89772E 

03 

-3.56484E 

03 

-7. 258558 

03 

4 * 27 52.7 E 

03 







-9,215006 

C3 

2. 157586 

03 

-3 .93 008 E 

03 

6.347 136 

03 

— 1 . 6962 7E 

03 

4.60163E 

03 

1*216 4 06 

05 

-7 . 16 539 E 

04 







-2 .353106 

04 

— 2*373616 

03 

4*647246 

03 

— 6 *53044 E 

03 

1.64 149E 

03 

L. 351786 

04 

-1.65569E 

05 

9 • 1 128 1 E 

04 







-6.469446 

03 

-4.880476 

01 

-3.711416 

02 

1*3729 56 

03 

5 . 982 70c 

02 

9. 24100C 

03 

-7. I6650E 

C3 

C.C0000E- 

-01 







-3.081326 

04 

- 1 • 1 1 456E 

C3 

6*035516 

03 

— 1*922 10E 

03 

5.14018E 

03 

8.009946 

03 

-1.531C6E 

C5 

6.63463E 

04 










MATRIX r C 2 ’ 8x8 


1*644486 

C3 

1.022486 

03 

6.57634F 

01 

— 4 * 8 3633 E 

01 

-1.677076 

02 

— 1 .552066 

03 

-2.894C7E 

03 

6. 54C9BE 

02 








03 

9 *00 3396 

02 

8. 103B3E 

02 

4. 37531 F. 

02 

1 • 444646 

02 

-3.267906 

01 

-1.119586 

-1.836126 

03 

4.473596 

02 









3 • 1 i 56 IE 

02 

5.20U5E 

02 

8.396476 

02 

3*330126 

02 

1.529326 

03 

—9 *911026 

02 

—1*756896 

03 

5.349776 

02 







-5 • 65560E 

02 

— 6* £5 3836 

01 

l .020496 

0? 

3.486726 

0? 

6.0C313E 

02 

— 2*5 1956 E 

02 

C* 137C3E 

C2 

3.805236 

02 









5.515386 

C2 

1 . 82 33 9E 

02 

5.23090E 

02 

-6.07225E 

02 

8.87192E 

03 

-6.498376 

03 

— { • 78829E 

03 

6.88676E 

02 




. 





-8.40573E 

02 

-5.34677E 

0? 

1*803016' 

02 

-1.37225E 

02 

-6. 32260E 

03 

1 • 3691 3 E 

04 

1.072C26 

C4 

—1.1391 86 

03 








03 

3.045836 

02 

1.854536 

02 

— 3 . 046 4 8 C 

02 

-3.21030E 

02 

-1.862986 

03 

2 • 51455E 

3*330136 

C3 

0.000006- 

-01 









—1*425416 

03 

— 1 • 1 2432 E 

03 

— 3 *412116 

02 

3.544166 

02 

—3*752076 

03 

5.72328E 

03 

7.716246 

03 

1.44666E 

04 
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MACH 1.2 


MATRIX , C 3 * 8 x 8 

1*1788 3 E 01 8 *202638—01 4.40056E 01 

1.03256F 01 -3.52417E 02 

4* 25572E CO 1.04562E 00 2*C7921E 01 

5.67483F 00 -2.02288E 02 

-5.2326CE OC -1 .4C556E 00 2.00159E 00 

3.12645F 01 -1.48572E 02 

-4.66615F 00 — 5*034 38F 00 -5.45414E 00 

—6* 95494 E 00 1.16358E 02 

— 1 * 79 866E 01 -6* 3092 3 E 00 -1.19362E 01 

2.C6944E C2 -1.01902E 03 

2 • 2 275 4 F Cl 1*91898E 01 2*930686 01 

— 1*9381 3F C2 1 • 25222E 03 

8* 10466E 00 3 * 37 7C9F 00 1.65893E 01 

-5.21456F C2 O.COOCCF.-Ol 
4*31 138 E-03 — 1 *34 863E 00 -1.10585E 01 
2.7C5C9E CO 1*30105E 03 


1*064C1E 01 2.18312E 02 -2*7i655E 02 

9.77275E 00 1.09550E 02 -1.49634E 02 

7.71825E 00 -3*33349E 01 -6.C7876E Cl 
-1.71370E 00 -4 « 54686E 01 5.49588E 01 

1 * 1 0855E 00 —2*6153 IE 02 -3.460C7E 01 
1.640C2E 01 2.87250E 02 -2*28288E 02 

-5 *66 132E 00 1.31 735E 02 -6.31852E 01 

-8*2 0690E 00 -2.26626E 01 I.C3947E 02 


MATRIX 'D*^ 8x8 

—6*01 543E— 01 1.25C28E 00 

-4.99637E 00 -2.02172E 00 
5 * 5£ 593E— 02 2.76872E- 01 

-2.29530E-02 -i.i9718E 00 
1.14 6 02F— 0 1 2 * C6 1 67F — 0 1 

— 2*23379 E 00 -1.03559E 00 

1* 5694CE-0 1 -2.25980E-01 
1.CE394E 00 4*900116-01 

-7.76629E-01 2.26827F. 00 

-1.66C68E Ci -6.2764BE 00 

— 1, 39204E-01 -1.48819E 00 

1.15399E 01 7.95609E 00 

-1.189156-01 -1.901616-02 
-7.89703F-C1 C.OOOOOE-Ol 

— 3*7359 IE-0 1 -1.183346 00 
-6.53250E CO 1.CQ902E 01 


MATRIX 'd x ' 1x8 

C^.C0C00F-03 3. OOOCOE-Q3 

3 . C0000E-03 3 . 00000E-03 


4.23218E 00 6.92064E-01 

5*4 1037E— 0 1 —1*3782 1E-0 1 
3 • 00507F— 0 l 3. 31097 E- 02 
0.76395E-O1 — fc * 3 8037 E- 02 
-7.76216E 00 1.74479E OC 

4.35370E 00 -7.96204E-01 
-2. 03780E-01 1*2 4083E— 02 

4.12853E 00 -fc. 36521E-01 


3.00000E-03 3.0C0C0E-03 


— 2 * 83688 E 01 -1.59500F 01 

— 5*4-5951 E 00 -1.22214E OC 
-1.18254F OC 1.52362E 00 

5.59713E 00 3.61562F OC 

—4 *24894 E 01 -2.57048E 01 
2* 52071E 01 i *422 246 01 

— 8 • 84319E— 01 — 1 • 97269 E 00 

2.BC687E 01 1.17362E 01 

3 • OOOOOE— 0 3 3 • C00C0E-C3 
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MACH X.* 


MATRIX *D 2 * 8x8 


-4.E7614E 

CO 

-3.471246 

00 

1.727276 

02 

3.23918F 

01 

-3.56392E 

00 

-1.31C93E 

00 

3.302206 

01 

1.91 31 96 

01 

-1.794366 

00 

-1.82138E 

00 

7.31272 E 

01 

1.657846 

01 

2.82827E- 

01 

6.505586- 

01 

—2.896366 

01 

-7.687636 

00 

-5.334166 

CO 

-1.083336 

01 

4 .953206 

C2 

5.58726E 

01 

1.5 28 ICE 

01 

I.2C789E 

Oi 

-4.112116 

02 

-1 .26776E 

02 

.1.056836 

00 

4.43741 E- 

01 

1.557CIF 

01 . 

O.OOOOOE- 

*01 

2.245196 

01 

6.57872E 

00 

— 5 • 7C £ 8 6E 

Cl 

-1.700946 

02 

MATRIX 

•d 2 

*1x8 


9.000006- 

03 

9.00000E- 

■03 

5.000006- 

■03 

9.00000E- 

*03 


-1.474816 00 8.94350E 00 

-3 •97599E 00 5.60154E 00 

-1.99934E 00 2.63242E OC 

-U37524E 00 -1.395C5E OC 
1.02551E 01 9.15584E 00 

1.42953E OI -8.16487E 00 
4.67606E-01 -7.65747E-Q1 
1*709 14E— 0 1 -7 *53121 E 00 


9. OQQOQF— 0 3 9.00000E-03 


-4.27613F 01 — 4 • EQ729E Cl 
-1.92228E 01 -5 • 2979 IE Oi 
-1.879616 01 — 1 • 556436 00 
2 • 6 8636E 00 -2.14256E CO 
— 5 • 17529 F 01 1.31025E 02 

8.46021 E 00 — 1.31313E 02 
-1.332066 01 1.1C278E 00 

2 • 6 1 C84E 01 —1*213966 Cl 


9 • 000006-0 3 9 • COOCOF-03 


MATRIX ’D^ 8x8 


1.5386C6 01 2.251976 01 

-6* 1 4 C i 66 C? -1.5318SF 02 
1.2755CE Cl 9.22*366^ 00 
-1.691466 C? -9.C5120E 01 
5.65656F 00 7.42 605E 00 

-3.48C47E 02 -7.851036 01 
1 * 6 C t C9E 00 -3.34C526 00 
7.55573E Cl 3.60098E 01 
3*9 171 2E 01 5.23573E 01 

-2.1U05E 03 -4.71510E 02 
-6.233516 01 -6.C7E346 01 
2.07157F 03 5.585206 02 

-5.5621CE 00 -i. 536036 00 
-6.33257E Cl O.OOOOOE-Ol 
-8.316UE Cl -3.9V591F Oi 
7.67320E 02 7.943716 02 


3.12849F 01 —3 *8 3271 E 01 
4.226206 00 -2. 69256E 01 
1.260386 01 -2.13696E 01 
1.277826 01 3.44589E 00 
7.239606 01 -6.524C2E 01 
7.967316 Oi 7.26057E 01 
? . 446446-01 3.76172E 00 
2.552G8E 01 5.22903E 01 


9.15851E 01 4.34563E C2 
1.57161E 01 2.74316E 02 
l • 0444 1 E 02 1.788136 02 
2.68556E 01 -7.04932E 00 
8.62461 E 01 3.331106 02 
3.961726 01 -1.510636 02 
4.36328E 01 -6.315416 01 
•1.403276 01 -3.12954E C2 


MATRIX *d^* 1- x 8 

l. 5C0OC6-O2 1 • 500C06— 02 l. 50000E-02 1.500CQE-02 1.500006-02 1.500C0E-02 


1 . 50C00E— 02 1 . 500CC6-02 
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MACH 1.2 


MATRIX *D4* 8x8 

- 1 . 87 AA IE 01 -2.8982AE 01 A.70626E 01 

0.2439CE 02 1.92697E 02 
-1.51C2EE 01 -1.391A6E 01 -2 . 6 8585E-0 l 
1.09868C 02 1.13859E 02 

-1.A7287E Cl -1.2CC29E 01 -8.5234AE 00 

1.55973E 02 9.87056E 01 

W*.C62?<iE 00 1*72 99 0E 00 -1.52033E 01 

-t*725A6F Cl -A * 55220E 01 

-5.A5781E 01 -7.Q2356E 01 1.AL562E 02 

I * 9 9 6 9 8 F C3 5.9384AE 02 

0. 629F3F 01 8.73C69E 01 -1.A9A66E 02 

-2.10603E 03 -7.54455E 0? 

1. A7298E Ci 4.A3855E 00 -2. 006776-01 
1 . 1 A 03 IE C2 0 • OOOQOE—O l 

9.6C750E Ci 5 • 69 8 1 8E 01 -5.81238E 01 
-1.00173E 03 -1.C0819E 03 


MATRIX *d^ r 1x8 

2 • ioOOGF— 02 2 • 10 OOOE— 02 2.100C0E-02 

2.100C0E-C2 2. 10000E-02 


3 *260716 01 2 *54066 E 02 -4*26439E C2 

2 *876506 01 8.21A1AE 01 -3.01765E 02 

2 #62 ? AOfc 01 — 3* 4 1 777 E 01 -2.79057F 02 
2 • A379CE** 01 -2*A8373E 01 -A*22A1CE Ci 
6*1851 3E 01 A • 04977 E 02 -5.44191E 02 
-8.716C3E 0l--5.73t.37E 02 5.25359E 02 

2*68570 E- 01 -I.37137E 02 5.24299E 01 

6 *089996 01 ~i*93887E 02 4*20399E 02 


2 * 1 OOOOE— 02 2.10000F-02 2.IOOOOF-0? 




3.0 SST WING MODEL CONTROL SURFACE MECHANIZATION 

This section discusses the analyses, design, and laboratory testing 
accomplished on the mechanization of control surfaces for the SST wing model 
flutter suppression system. Actuation system components have been identified 
to provide the desired performance, but work remains in incorporating the systems 
into the model. 

3*1 Introduction 

Since 1970, Boeing -Wichita has provided support to the NASA- Langley 
Research Center Aeroelasticity Branch on a research program to demonstrate in the 
Langley transonic dynamics wind tunnel an active flutter suppression system on a 
l/l7 scale SST semispan .wing model. The initial effort was directed toward mech- 
anizing an electromechanical actuation system for the model trailing edge control 
surface (see Reference 2). This system was subsequently determined inadequate 
for the flutter suppression system demonstration, and in 1971 an effort began to 
develop electrohydraulic actuation systems for the model midspan leading and 
trailing edge control surfaces. The following paragraphs describe the work that 
has been accomplished on this development. 

The design philosophy followed in this development is straightforward: 
establish performance requirements that will lead to a successful wind tunnel de- 
monstration of the flutter suppression system, purchase and/or design and fabricate 
components that satisfy the performance requirements, and breadboard test the 
systems thoroughly before installation in the model. 

The NASA- developed flutter suppression system was formulated on the 
assumption that ideal actuators would be used to drive the control surfaces. 

Based on this assumption, a dynamic requirement of minimum phase and gain variations 
was set in the frequency range 5 to 25 Hz, with the model flutter mode frequency 
about 12 Hz. A goal of no more than 15 degrees phase lag at 25 Hz was chosen. 

TTie actuation systems must be capable of ± 10 degree amplitude up to 
25 Hz, and produce at least 20 inch-pound maximum torque. No external leakage of 
hydraulic fluid can be permitted in either the model or the wind tunnel test sec- 
tion. The hydraulic actuators must mount close to the control surface hinge lines 
without violating the wing airfoil. 

A subminiature rotary actuator was designed and fabricated to be used 
with high performance Moog servovalves to actuate the midspan leading and trailing 
edge control surfaces. Detailed discussions of the design, testing, and analyses 
accomplished to date are presented in Sections 3.2 and 3.3. Section 3.4 presents 
a discussion of work necessary to couple te the control surface mechanization. 

Section 3.5 contains nomenclature and detail drawings. 
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3*2 Control Surface Actuation System Design 

A survey of current "off the shelf" subminiature eleetrohydraulic 
components indicated that, while satisfactory servovalves were available, special 
actuators would have to be designed and fabricated to meet the requirements. 
Therefore, two single vane rotary actuators with 0.5 inch vane radius and length 
were designed and fabricated for installation at the control surface hinge lines 
without protruding into the air stream. 

3.2.1 Actuator Design 

The design torque requirement of at least 20 in- lb maximum was based 
on control surface hinge moment estimates. To insure satisfactory operation of 
the actuation systems during the wind tunnel tests, a design goal of achieving a 
higher maximum torque capability was set, within the geometric constraints of the 
model. A limit of 1000 psi was chosen for hydraulic supply pressure to minimize 
fluid sealing difficulty in the actuator and lines. 

A detailed drawing of the final actuator design is shown in Section 
3.5. The actuator body and end caps are 2024 aluminum, and the vane and shaft 
are 17-4PH high strength stainless steel. Vane sealing is provided by coating 
the vane with adiprene in conjunction with a straight segment of an AN 5227 
0-ring behind the shaft. Gasket seal between the end caps and actuator body is 
accomplished with two loops of 3 mil copper wire set in a light adhesive arcund 
the vane cavity. The l/8-inch precision shaft is supported by ball bearings in 
each end cap and AN 5227-1 0-rings provide fluid seal around the shaft. A photo- 
graph of one of the actuators prior to assembly is shown in Figure 3*1. 

The two actuators were designed to be identical, with the actuator 
body and end cap thickness machined to lie within the airfoil at the inboard edge 
of the leading edge control surface. Layouts of a proposed installation at each 
surface are shown in Section 3.5. Special l/8-inch tube fittings were designed 
and fabricated to provide 0-ring seals at the actuator ports. 

The actuators were designed to produce 4l.O in- lb torque with two- 
thirds of the 1000 psi supply pressure across the actuator vane. Usual actuator 
design procedures were followed (see Reference 3)- A method of attaching the 
surfaces to the actuator shaft has not been determined. The use of a tapered 
pin would reduce the stress area of the shaft and aluminum tubing the surfaces 
are mounted on to the point where failure would occur at 3O-35 in- lb. 

3.2.2 Eleetrohydraulic Servovalve Selection 

Servovalves were selected to produce 10 degree actuator amplitude 
capability up to 25 Hz in a no-load condition. The servovalves selected are 
Moog Series 30 flow control type (part number 030 A 17010 E 022 F4) rated at 
1.7 in^/sec no-load flow rate at 1000 psi supply pressure. With a load pressure 
of two- thirds the supply pressure, the maximum flow rate is 58$ of the no-load 
flow rate. With the design torque load (41.0 in- lb) on the actuator, the servo- 
valve can produce 10 degree actuator amplitude up to 14 Hz and 5.8 degrees at 
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25 Hz. This capability was estimated with the assumption that the servovalve and 
actuator would be mounted close together, whereas in the model they will be sep- 
arated by about 42 inches of stainless steel hydraulic tubing. Fluid compliance 
will degrade this estimated performance slightly, but the actual performance 
should be more than adequate to work the 12 Hz flutter mode. 

The Moog Series 30 servovalves were selected primarily because of 
their high bandpass dynamic capability. The apparent undamped natural frequency 
(-90° phase point) is about 240 Hz with 0.50 damping ratio. 

Port manifolds, we.re designed and fabricated from 2024 aluminum to use 
with the servovalves (see Section 3.5)» 0-ring seals are used between the mani- 
fold and the servovalve ports. The manifolds (and servovalves) will mount on the 
wing mount plate within the fuselage fairing. 

During the wind tunnel tests, the flutter suppression system will be 
mechanized on an analog computer. The error voltage to drive the servovalves 
(and subsequently the control surfaces) will be formed as the output of a computer 
operational amplifier. A servo amplifier has been designed and built using a 74l 
operational amplifier to accept the error voltage and produce current to drive the 
servovalve as a constant function of the error voltage, independent of frequency. 

A schematic of the servo amplifiers is shown in Figure 3*2. The servo amplifier 
will permit maximum current to the 500 ohm servovalve coils (in series) up to 
about 25 Hz. Thus, the maximum flow rate capability of the servovalves can be 
attained up to this frequency* The two servo amplifiers are mounted on a small 
circuit card and will mount close to the servovalves under the fuselage fairing. 

Figure 3.3 shows a simplified plumbing diagram from the hydraulic 
supply pump to the servovalves, and return. The pump must supply MIL-H-5606 
petroleum base fluid at 1000 ± 50 psi at a maximum flow rate of 3*4 inVsec (.88 
QPM). A 5 micron nominal, 15 micron absolute filter is required in the pressure 
line. A 25 in3 (or larger) accumulator, pre charged to 5 00 psi, in the pressure 
line is recommended to damp out pressure transients or ripple from the pump , and 
to keep the pump from having to respond to high load frequencies. 

3.2.3 Actuation System Installations 

A sketch of the actuation system installations is shown in Figure 3.4. 
Two feet of l/8-inch outside diameter stainless steel tubing will be used in each 
line, from the actuators inboard. The remaining distance will be l/4-inch tubing. 
This choice of tubing sizes represents a compromise between pressure loss due to 
fluid friction at maximum flow and compliance in the lines between the servovalves 
and actuators. The hydraulic tubing will be set into- channels cut into the balsa- 
wood, and bonded with epoxy to the aluminum alloy plate out to about one foot from 
the actuators, where an elastic adhesive will be used. The channels will be filled 
in over the lines before the wind tunnel tests. 


Layouts showing details of the actuator installations at the leading 
and trailing edge control surfaces are included in Section 3*5* The actuators 
will mount on the model aluminum alloy plate. The actuator shafts will be attached 
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directly to the aluminum tubing the balsawood control surfaces are mounted on, 
and will provide the inboard support for the surfaces. The outboard supports 

shown in the layouts will be used for laboratory testing of the actuation systems, 

but during wind tunnel tests pivot bearings must be used to prevent binding as 
the wing flexes. 

An angular position sensor is being designed to mount on the actuator 
shaft on the side opposite the control surface. The sensor will use two silicon 
photocells mounted on a common brass base, with 0.010 inch gap between the cells. 
A semicircular area of light will be projected on the cells such that the differ- 
ential voltage generated is proportional to the shaft angular position (see 
Reference 2). 

3*3 Baseline System 

3.3.1 Baseline System Testing 

A baseline actuation system was assembled for breadboard testing to 
evaluate the hydraulic actuation systems. The baseline system, shown in Figure 
3.5, has about 15 inches of l/4-inch outside diameter stainless steel tubing and 
4 inches of l/8-inch tubing between the servovalve and actuator. A d.c. film 
type potentiometer provides the position feedback signal and the feedback loop 
is closed on an EAI TR-48 analog computer. Commands from the analog computer 
pass through the servo amplifier to the servovalve. The trailing edge control 
surface is used in the baseline system because it possesses higher inertia than 
the leading edge surface, resulting in a lower frequency surface-shaft mode. An 
MS 28797-1 accumulator is used in the pressure line to isolate the system from 
the pump in the hydraulic test bench. 

The baseline system with position feedback only is unstable with the 
control surface connected to the actuator shaft. The instability appeared during 
testing as a sustained 55-9 Hz, 1.2 degrees peak to peak oscillation as the supply 
pressure was increased from zero at about 550 psi, with a nominal loop gain of 
633«7/sec. This gain had been determined through testing of the baseline system 
without the control surface to provide less than 15 degrees phase lag at 25 Kz 
(see frequency response shown in Figure 3*6). The frequency response indicates 
a lightly darped peak at about 99 Hz. No attempt was made to determine the peak 
exactly during the laboratory testing due to the high auplitudes near the resonant 
frequency. 


3*3,2 Baseline System Analysis 

A sinplified, linear mathematical model of the baseline actuation 
system was developed using the system sketch shown in Figure 3.7 and the frequency 
response of Figure 3*6. This mathematical model will be used to predict additional 
feedback compensation required to stabilize the baseline system. 

Consider the sketch shown in Figure 3.7. This combination free body - 
flow diagram is drawn to represent the baseline system. Summing torque about the 
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actuator shaft produces the equation 

Ie * + D "» 3F + k s®a- k 3 q s = t a 

where Ip, Q is the equivalent rotary inertia of the hydraulic fluid and the actuator 
vane, Djjq is an equivalent linear (viscous) danping coefficient and T^ is the torque 
developed by the actuator. The developed torque can be determined by assuming the 
pressure acts on the vane area at the average radius: 


Ta = (P, - P 2 )[(v„- Vs)*] = i-(f? - P a \(v- D 2 - r*)jL = C a {P,-F^ 

2. c ~ 

where C A is the actuator effectiveness coefficient with units of in^. 

The Moog servovalve dynamic characteristics can be assumed constant 
for low input commands and loads. The approximate servovalve equation of motion 
can be expressed as 


<H? 


+ 2-t v U) n 


dQ(4) 

v d-t 


+ u>* v Qtt> 


Ky Lit) 


where i(t) is the coil current and Ky = ^ na3C /^may ^ jOW ' rA ^ e hydraulic 

fluid to the servovalve divides into dating flow (Qj)) due to leakage across the 
spool, flow due to fluid compressibility (Qh)> suid actuator flow (Q^). The damping 
flow is accounted for in the approximate servovalve equation given above. The flow 
in the line from the servovalve to the actuator is 


Qti) = <?««.) + q a w . 


The compressibility flow can be determined from the definition of the 
fluid bulk modulus 

AP $ 

& = av7v or V = AP 

where AV is an incremental change in fluid volume due to an incremental change in 
pressure, AP, and V is the fluid volume on either side of the actuator vane. Both 
sides of the equation can be divided by an incremental time, At, and in the limit 


tim f_£ AVT _ lirr, AP 

At-*0 L V A-t J “ at-*o At 

produces 

A 

v a-t ~ d-t n or at 


£ 

v 


Quit) . 
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The actuator flow can "be obtained by recognizing the flow rate is the 
rate of change of volume in the actuator on one side of the vane due to an angular 
displacement, 9^. 



Then, 

Q(t) 


- dP 

& d* 


+ C, 


<AQa 

d*t 





The load pressure P 3 .-P 2 I s twice the pressure AP since pressure increases on one 
side of the actuator vane by AP and decreases by the same amount on the other side 
as the servovalve opens. Then the flow equation may be written as 


Q 


_V_ dt'P.-Pfc) . r d©*, 
“ 20 Si + IT 


This equation can now be used to obtain the actuator torque 
tion of servovalve flow rate. Integrating from some time t Q to time t 
form 



2 0 


C A ©A 


as a fUnc- 
gives the 


where it is assumed (P^-P 2 )(t 0 ) = 0 and 6 a(^o) = 0* From the torque equation above, 

t a c« - & c A fQ(r)dr - ^ cl © A(fc ) . 

**0 

Substitution of this equation into the actuator equation of motion gives 


or 


- d Z ®A 

» 

t d" e A 


+ D de* , 
D *» At 


*9 


<H- Z 


« s e A -K 3 0 s = y C A fotnclT - ^ c* e, 
+ D«» ^ + (v c a + K 3 )e A - K s ©s = fc A Jfl(T,dr. 

t* 


The other equation of motion required to describe the baseline system 
dynamic behavior can be obtained by summing torque about the surface shaft. 


- Ks ®a + *s ^ + ®s « o . 

Assuming zero initial conditions, the Laplace transformation of the two equations 
produces 
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(t w S 2 + 0 £q S + Y c| 4- K s \e A <3) - K s 0s ts) = ^ C A (-g-) Q(S) 
-« S © A (S) + (l s 3 2 4 K^e^S) = O 

The transfer function 6 a (S)/q(s) can be determined by applying Cramer’s rule 

« egc^s* + k 5 ) 


<S) = 


3 [(!„<? 4 D w s + & C^KsVl s S 2 <■ Ks'j - Kj] 


This transfer function for the baseline system without the surface 
inertia can be obtained by setting Ig — 0 and Kg — 0 in the transformed equations, 


^ (S') = 


gc A 


s[. 1 eoS Z -+ Deq S +• %§■ c| J 


The block diagram shown in Figure 3.8 represents the baseline actuation 
system without the control surface. The frequency response of this system, Figure 
3*6, indicates a resonant peak at about 99 Hz and 0*30 damping ratio. This mode is l 
the closed loop equivalent hydraulic fluid-actuator inertia mode. The closed loop 
transfer function is 

©A (s . = 57.3K AMP K v tog v (2 6C A /VXg^ ( DEG /VOLT) 

Vc s(S 2 + 2J,U>„ V S + <)(S 2 + |S3 S + 1|^} + 57.3^K AMP K v w^(lg) 

where the fifth order denominator has two pair of complex conjugate roots and one 
real root. The closed loop denominator can be written in terms of open loop 
quantities: 

(s z 4- 2 T, *>„, 3 + «*>«,)( S 3 + q, S l 4 q*S 4 q s ) = s(s l + 2 * v u>r v S + <*>£)(s* + gj S + 

’ + S7.3K P K M ,p)<vW^(l^fi) 

where ^ .01 and W ri| « (2tf)(99). Five linear simultaneous equations in the 

five unknown coefficients (o, » , and — 1 can obtained 

3 XgQ '•*tQ * . 

by substituting the known quantities in the above equation and equating coefficients 
of like powers of S. The solution of these equations gives the open loop hydraulic 
flu id- actuator inertia mode as 

s 2 * + 2*9 3 4 - 2 ^ Ca 3 S e -+- 4 . 7 TT Xl 0 3 s ■+ 3.080 X lO 6 * (S+ 768 . 4 XS+ 4 O 0 S. 3 ) 
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FIGURE 3.9 

BLOCK DIAGRANA OF BASELlME SYSTEM WITH CONTROL SURFACE 









indicating that this mode is overdyed. This appears reasonable due to the 
friction between the actuator vane and body, and between the shaft and O-rings. 

The equivalent danping and inertia are ^ = 4.012 and l EQ = .0008^ in-lb-sec . 

The transfer function 6 ^(s)/q(S) with the control surface included 
(derived above) may be written in the form 


£* (s) = 


USt (S* + ~) (57.3) 

VIu 


/ DEG \ 
\,.wVSECl 


s f / s 2 . Oso s + 2gC *- * — \(s 2 + -)- 1 

*- ( iTo vl eQ Is 


The shaft spring rate was estimated at 368.9 in- lb/rad, considering noth the actua- 
tor shaft and the aluminum tubing that the surface is mounted on. The hydraulic 
fluid- actuator inertia mode becomes 

ft «?££?• K« _2 - Jb ^ & 3G8.9 IN-LB/RAO 

5 l + ^ s + Ipft- - ~ = S + 4.T77* 1 0 3 S +3.030x10 + 


ra 




= S 2 + 4.777 Xl0?5 + 3.519 Xl0 6 
= (5 4*9t0.l)(S + 386 G.G) . 

The midspan trailing edge surface inertia was supplied by NASA and is 9*^32 x 10 
in-lb-sec^. The surface inertia- shaft mode is 


s 2 + hs = 


S 2 4 


3&8.<* H4-t©/RA0 


= 5 z + 3.qtl XtO 5 * 


I 5 * <h432 * 10“ 4 in-i.s-sec z 

The transfer function denominator then becomes 

(3 €>S . 9 ) 

s[.( 5 2 +4.777X10*5 +3.5>I9x|0 6 )(s 2 +3.911 XIO 5 ) - (ooog4)( 


I A \ f rtATlOAO) 


-] 


- 5^s 4 + 4-TT7 S 3 + 3.910 xl0 6 S 2 + I.®&&X10 9 S + 1.204 xIO 12 ) 
= S(S -(- 859. 3)(s + 3870.8)(s+ 23.43 tj601.a) 


and the numerator is 


ZB C A 

VX E Q 


(57.3)(S Z + j?) 


20.9 K \Q?L&(tN*-)(*0Gl5 3 0f<?/R4p) 3 qn XfO ff ) 
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2.870 xIO 9 (S 2 4 3.9 If *ID 5 ) . 
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The block diagram of the baseline system with the control surface 
included is shown in Figure 3*9, with the position feedback loop closed. Addi- 
tional feedback condensation must be determined to stabilize the system, and the 
analysis extended to the actuation systems as installed in the model. The com- 
pensation that stabilizes the baseline system should also stabilize the actual 
systems. The compensation will be incorporated into the breadboard system and 
performance verified through laboratory testing. 

Several simplifying assumptions were made in the above derivation. 
The most significant is the assumption that the servovalve behavior can be des- 
cribed by a constant coefficient, linear ordinary differential equation. The 
change in flow gain due to load pressure has been neglected, but since this is 
reflected as a loss in loop gain, the mathematical model is conservative. The 
servovalve apparent undamped natural frequency and damping ratio were taken from 
a Moog brochure for the Series 30 servovalve. Combining the hydraulic fluid mass 
with the actuator vane mass into an equivalent inertia is not new (see Reference 
4). The correctness of all assumptions will be determined through laboratory 
testing of the breadboard system, with the analytically identified compensation 
incorporated. 

3*4 Remaining Work 

Components have been selected for the electrohydraulic actuation 
systems, but work remains in integrating these components into stable systems 
meeting all the performance requirements. The primary item to be accomplished 
is the determination of feedback compensation (in addition to actuator shaft 
angular position feedback) required to stabilize the systems. A mathematical 
model of the baseline system has been developed to permit analytical determina- 
tion of the required compensation. The compensation will be tested on the bread- 
board baseline system to verify the analysis results. The breadboard system will 
be modified to include the line length dictated by the wing geometry. When the 
systems have been proven through breadboard testing, they will be installed in 
the model and retested to verify performance. 

Some effort must also be spent in perfecting the silicon photocell 
angular position sensor. The concept has been proven through laboratory testing, 
but the design for this particular application must be perfected to provide a 
reliable actuator shaft position feedback signal. 

Provisions have already been made to install the midspan trailing 
edge surface actuator in the model, but nothing has been done to the model toward 
installing the leading edge surface actuator. 

3»5 Nomenclature and Drawings 

This section includes nomenclature and detail drawings of fabricated 
hardware required for the control surface actuation systems. 
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3.5.1 


Nomenclature 


Symbols used in the mathematical development in Section 3.3*2 are 
listed below. 


STOOL 

DEFINITION 


UNITS 

C A 

Actuator torque effectiveness coefficient 
(C A = 1/2 (r^-r^JA) 


m3 


Equivalent viscous dancing coefficient of actuator 
and hydraulic fluid 


in-os 

rad/sec 

i 

Servovalve coil current 


amp 

X EQ 

Equivalent rotary inertia of actuator vane and 
hydraulic fluid 


in- lb- sec 2 

is 

Rotary Inertia of control surface with respect to 
its hinge line 


in- lb- sec 2 

Kp 

Actuator shaft angular position feedback gain 


volt /deg 

KS 

Torsional spring constant of actuator shaft and 
surface tubing 


in- lb /rad 

K V 

No-load flow gain of servovalve 


in3/sec 

anp 

JL 

Length of actuator vane (= .50 inch) 


inches 

p 

Hydraulic fluid pressure 


1'o/in 2 

Q 

Hydraulic fluid flow rate from servovalve 


in^/sec 

r D 

Actuator vane radius, measured from shaft axis 
(= .50 inch) 


inches 

r S 

Actuator shaft radius 


inches 

S 

Laplace transform operator 


l/sec 

Ta 

Torque developed by actuator 


in- lb 

V 

Volume of hydraulic fluid on one side of actuator 
vane (from servovalve to vane) 


in3 

p 

Bulk modulus of hydraulic fluid (taken as 
1.9 x 10 5 lb/in 2 for MIL-F -5606 fluid) 


lb/in 2 
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SYMBOL 


e A 
©S 
«j n v 


DEFINITION 

Apparent damping ratio of open loop servovalve 

Actuator shaft angular displacement 

Control surface angular displacement 

Apparent undamped natural frequency of open loop 
servovalve 


3.5.2 


Detail Drawings 


UNITS W 


radian 

radian 

rad/sec 


Copies of detailed drawings of the actuator (2 sheets) and servovalve 
port blocks, and layouts of the actuator installations in the model, are presented 
in this section. 
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4.0 B-52 AEROELASTIC MODEL RIDE CONTROL SYSTEM SYNTHESIS 

The synthesis of a ride control system for the NASA one-thirtieth 
scale B-52E aeroe las tic model is described in this section. The final system 
produces more than 30 percent reduction in rms vertical acceleration (airplane 
response to random atmospheric turbulence) along the entire model fuselage, 
using forward body horizontal canards, inboard wing flaperons, and elevator 
control surfaces. This system will be mechanized on the model and demonstrated 
in the Langley transonic dynamic tunnel, 

4.1 Introduction 

A three phase study to design a ride control system for the B-52 
aeroelastic model was formulated in 1970 and reported in Boeing Document D3- 8390-2 
(Reference 5). This document is included as Section 3 of D3- 8390-4 (Reference 2). 

The objective of this study was to design and evaluate a ride control 
system (RCS) for the B-52 aeroelastic model with maximum performance for minimum 
model modifications. During Phase I, the ride improvement attainable using the 
existing elevator and aileron control surfaces was evaluated. These surfaces 
were capable of producing no more than a five percent reduction in rms vertical 
acceleration (a) at the pilot station (BS 172) for a random vertical gust dis- 
turbance. Performance attainable with canard and canard/elevator systems was 
investigated during Phase II, Satisfactory A reduction was obtained at the pilot 
station, but mid. and aft body reductions did not meet the design goal of thirty 
percent. A ride control system was designed and evaluated during Phase_III using 
canard, flaperon, and elevator control surfaces. This system produces A reductions 
in excess of thirty percent over the length of the model fuselage. 

The remaining work includes using updated cable-mounted model equations 
of motion for evaluating the system performance, refining the system design and 
for generating system responses for correlation with wind tunnel data. 

4.2 Design of Ride Control System 

The objective of this study was to analytically demonstrate feasibility 
of improving passenger/pilot ride using active controls. The control surfaces, and 
surface locations, considered in the study are shown in Figure 4.1. The system 
design goal was to obtain a minimum of 30 percent reduction in fuselage rms verti- 
cal accelerations In a random atmospheric turbulence environment. 

The design study was conducted using model equations of motion that 
did not include the cable mount effects on the model dynamic behavior. All data 
shown in the following pages is model scale, but wing and body station designations 
are airplane scale. The mathematical model included the first 14 symmetric degrees- 
of-freedom for wind tunnel test condition 1 (Reference 6). The model equations 
included lift growth effects and the Von Karman gust spectrum. The characteristic 
gust length used in Phase I and II was 16.67 feet and 83.3 feet (equivalent to 
500 feet and 2500 feet in airplane scale) was used in Phase III. Phase III was 
believed to be more representative of actual flight conditions. Actuator dynamics 
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B-52 control surfaces and locations 


figure 4.1 


in Phase X and II were represented by a first order lag, while Phase XIX included 
the actual second order actuator dynamic behavior. 

The RMS-FSD analyses conducted in this study included the full 0-40 
cps frequency range of the model equations. All percentage reductions in vertical 
acceleration presented in this section are based on the full frequency range. 

4.2.1 Phase I - Existing Control Surfaces 

The B-52 aeroelastic model was originally constructed with elevator 
and aileron control surface actuation systems. The objective of this phase was 
to determine the ride improvement attainable with these control surfaces. 

Figure 4.2 shows the effects of elevator and aileron ride control 
systems on rms vertical acceleration along the model fuselage. Corresponding 
results tabulated as percent change in rms vertical acceleration are given in 
Table I. The results show that the elevator can significantly reduce mid and 
aft body vertical accelerations, but the ailerons are considerably less effective. 
The two systems were synthesized independently to produce the maximum vertical 
acceleration reductions attainable with the elevator and ailerons. Block diagrams 
of the elevator and aileron ride control systems are shown in Figures 4.3 and 4.4. 
Results with both loops closed simultaneously show that the elevator system alone 
is better. Since the aileron system alone was not effective in reducing fuselage 
vertical acceleration, no attempt was made to refine the aileron system with the 
elevator loop closed. 


TABLE I 

EFFECTS OF ELEVATOR AND AILERON RIDE CONTROL SYSTEMS 
ON RMS VERTICAL ACCELERATION 




in RMS Vertical Acceleration I 

Body 

Station 

Elevator 

RCS 


Elev. & Aileron 
RCS 

172 

-4.4 


+3.3 


-1.0 

860 

- 16.7 


-1.7 


— 11 . 4 

1655 

— ■ 

- 20.7 


-9.3 


-20.2 


Figures 4.5 through 4.l6 show PSD and RMS data which support the 

results given in Figure 4.2 and Table I. Some of the accumulative RMS plots are 

not complete out to 40 cps, but under each plot a final RMS value (corresponding 

to 40 cps) is printed. The units for the PSD and RMS axes respectively are 

(g^/ft/secj^/rad and g f s/ft/see. These plots are all with respect to a 1 ft/sec 
RMS vertical gust on the model. Due to velocity scaling, a 1 ft/sec gust lias the 
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same effect on the model as a 5*48 ft /sec gust has on the full scale aiiplane. 
Equivalent airplane values can be obtained by dividing the model RMS values by 
5.48. Ride performance was evaluated at three stations along the fuselage: 

BS 172 (pilot’s station), BS 860 (approximate eg), and BS 1655 (aft fuselage). 

The PSD’s in Figures 4.5 through 4.8 show that the model gust response 
at BS 172 has contributions from rigid body motions and several structural modes 
and that neither the elevator nor the aileron ride control system changes the 
response significantly. 

Figures 4.9 through 4.l£ show that the model gust response at BS 86 0 
and BS 1655 are similar in that they both indicate predominant rigid body and 
first structural mode contributions. The elevator ride control system is effect- 
ive in reducing rigid body response but neither the elevator nor the aileron 
effectively couples with the first structural mode. 

Pitch rate feedback from BS 86 0 is used in the elevator ride control 
system as depicted in Figure 4.3. Its effect on model dynamics may be observed 
in the root locus of Figure 4.21a. Primary coupling is with the short period mode. 
Figures 4.17 and 4.18 indicate model elevator displacement and rate in a gust en- 
vironment. Equivalent airplane surface activity may be obtained by dividing dis- 
placement by 5.48 and rate by (5.48) . 

The aileron ride control system shown in Figure 4.4 uses acceleration 
feedback from the wing (WBL 487*5 which is near the aileron). Root loci for this 
system used alone and in support of the elevator system are shown in Figures 4.21b 
and 4.21c. Aileron surface activity is shown in Figures 4.19 and 4.20. 

The 12th structural mode showed no coupling with the systems and was 
omitted from the root loci in Figure 4.21. 

The systems presented have at least 6 db gain margins and 60 degrees 
phase margins. 

4.2.2 Phase II - Canard Surfaces 

Phase I results indicated that a control surface on the forward body 
would be necessary to improve ride at the pilot station. For preliminary analysis 
purposes, a point force (P.F.) representation of canards at BS 172 was added to the 
equations of motion. The canard coefficients were obtained by transposing the BS 
172 Z-modal coefficient row. 

Phase II included evaluation of a canard ride control system and a 
canard/elevator system. Effects of these ride control systems on RMS vertical 
acceleration are given in Table II . 
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TABLE II 


<r 


EFFECTS OF CANARD AND ELEVATOR RIDE CONTROL SYSTEMS 
ON RMS VERTICAL ACCELERATION 


Body 

Station 

Percent Change in RMS Vertical Acceleration"] 

Canard (P.F. ) 
RCS 

Canard (P.F.) and Elevator 
RCS 

172 

-37.9 

-35.0 

860 

-10.2 

-23.7 

1655 

+ 1.5 

-12.2 


7 , at BS 172 was sensed and fed back to the canard through the feedback 
transfer function 


-.0001s 

(•00637s + l)( *02S + I) 2 . 

The elevator system was the same as presented in Section 2,1 for Phase I. 

Figures 4.22 through 4.27 show the PSD and RMS data from which Table 
II was derived. The canard system provides adequate vertical acceleration reduc- 
tion at the pilot station and some reduction at the mid body station. Comparing 
Figures 4.14 and 4.27, it is evident that the canard excited the aft body 1st 
structural mode and the 10 cps mode. With both the canard and elevator loops 
closed, the reduction at the mid and aft body is lower than the reduction obtained 
with the elevator system only. The canard system was designed independent of the 
elevator to produce the maximum vertical acceleration reduction possible with the 
forward body canards. 

Figure 4.28a presents a root locus of the canard system and Figure 
4.28b shows the effect of the elevator system on the closed loop canard system. 

The canard system is a wide bandpass system that causes large movements of the 
higher frequency structural mode poles. 

None of the systems presented thus far significantly improved the 1st 
mode response at the mid and aft body stations nor indicated adequate potential 
for reducing vertical accelerations by 30 percent over the length of the fuselage. 
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Phase III - Canard- Flaperon-Eleva tor Surfaces 1 

The results of Phases I and II indicated that a three surface control 
system would be required to achieve the design goal of 30 percent vertical acceler- 
ation reduction over the length of the model fuselage. Table III presents the per- 
cent reduction obtained with a canard- flaper on-elevator system. 

TABLE III 


EFFECTS OF CANARD- HoAPERON- ELEVATOR RIDE CONTROL 
SYSTEM ON RMS VERTICAL ACCELERATION 


Body 

Station 

Percent Change in RMS 
Vertical Acceleration 

172 

-36.4 

510 

-35.6 

860 

-39.2 

1237 

-32.7 

1655 

-31-0 


These reductions were computed using the Von Karman gust spectrum with a character- 
istic length of 2500/3O = 83.3 ft* The numerical values as a function of body sta- 
tion are shown in Figure 4.29. Figures 4.30 and 4.34 present the RMS-PSD data 
supporting Table III. 

Figure 4.30 indicates that the short period mode has been effectively 
suppressed at the pilot station, the 1st structural mode shows an insignificant 
increase and the 2nd through 10th structural modes have been suppressed. 

Aft from the pilot station, the basic airplane short period and 1st 
structural mode response increase mono tonically and the higher frequency structural 
mode responses become nearly insignificant by conparison. Using the elevator and 
flap system, the short period and 1st mode can be controlled. The increased res- 
ponse shown by Figures 4.31 through 4.34 in the frequency range of 4 to 12 cps 
results from the flap and canard systems. 

For Phase HI, the mathematical representation of the canard was up- 
dated to include canard aerodynamics. Aerodynamic flaperon representation was 
also used. Second order actuator dynamics based on control surface moments of 
inertia were used as shown in Figure 4.35* A significant difference between this 
system and the previous systems is that the elevator is driven with Z sensed at 
the aft body station which permits significant reduction of the 1st structural 
mode response. 
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Washout filters were included in each of the feedback, systems to pro- 
vide satisfactory handling qualities. Normal acceleration and pitch rate due to 
step elevator commands are presented in Figures 4.36 and 4.37 for the model with 
and without the suppression system. With the suppression system operating, the 
model response is slightly degraded but still considered acceptable. 

The root locus of Figure 4.38 indicates the reason for the response 
degradation. The notch filter pole associated with the elevator loop moves closer 
to the origin than the short period mode, and thus becomes the dominant system 
root. The filter was introduced to provide lag at the short period frequency and 
lead at the 1st structural mode frequency (- and + 15 degrees respectively). Both 
the elevator and the flap loops are rolled off at relatively low frequencies. The 
canard loop, however, has a wide bandpass to make possible the vertical accelera- 
tion reduction for modes 2 through 10 shown in Figure 4.30* Additional compensation 
was introduced into the canard loop to provide the required gain margin of 6 db. 

The movement of the canard actuator root sets the upper limit for the canard loop 
gain. Nominal gains are those shown on the block diagram of Figure 4.35. 

RMS control surface displacement and rate requirements per ft /sec gust 
are presented in Table IV. Equivalent airplane scale values can be obtained by 
dividing RMS displacement by the velocity scale factor, and RMS rate by the velo- 
city and frequency scale factors. For example, in airplane scale, the canard 
rate would be l6l.8 - ^ dee/sec 

5 * 39 n7ifr . 


TABLE IV 

RMS CONTROL SURFACE REQUIREMENTS 


Surface 

Displacement 

deg 

Rate 
deg/ sec 

ft/sec 

ft/ sec 

Canard 

3-77 

161.8 

Flaperon 

2.16 

28.2 

Elevator 

0.34 

9.16 


Figures 4.39 through 4.4l present the RMS-PSD plots from which Table 
IV was derived. Figures 4.42 through 4.47 present Bode magnitude plots of the 
control surface displacements and rates for 1 ft /sec sinusoidal vertical gusts. 

For reference, the RMS values corresponding to Figures 4,30 through 
4.34 are tabulated as a function of frequency in Tables V. and VI. All percentage 
reductions quoted in this report have been based on the frequency range 0 to 40 
cps. 
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Figure 4.48 shows the effect of reducing all feedback gains to .7 of 
the values given on the block diagram (Figure 4.35). Note that the slope of 
the curve for the aft body is considerably less than for the other two locations. 

It is estimated that operating at .7 of the normalized gains would reduce canard 
and flap displacement requirements by 22 percent and elevator displacements by 
27 percent. Decreasing the gains to .7 increases the response of the short 
period and 1st structural mode as shown in Figure 4.49. At the same time, high 
frequency structural mode responses are increased on the forward body and decreased 
on the aft body. The dotted curve shows the effect of increasing the normalized 
elevator gain to 1. with the canard and flaperon gains remaining at .7. A summary 
of the percent RMS reductions associated with these gain variations is presented 
in Table VII. 


TABLE VII 

EFFECTS OF RCS FEEDBACK GAINS 
ON RMS VERTICAL ACCELERATION 


Body 

Station 

Percent Change in RMS Vertical Acceleration 
NORMALIZED FEEDBACK GAINS* 

Kq = 1 . 
Kp - 1 . 

k e = 1. 

Kc = .7 

% = .7 

K e = .7 

c — t — i—4 

« • 

it ti 11 

& j?jp 

172 

-36.4 

- 29*5 

-28.8 

860 

-39.2 

-32.3 

-33.8 

1655 

-31.0 

-28.0 

-31.8 


4.3 Remaining Work 

A ride control system has been synthesized to produce at least 30 
percent reduction in vertical acceleration along the entire fuselage. This 
system was designed using model equations of motion that did not include the 
cable mount effects. The cable-mounted model equations of motion will be up- 
dated to account for the change in cable attach point, and will then be used to 
evaluate and refine the RCS before the wind tunnel tests. 

A procedure must be established to bring the RCS on line on the model 
in the tunnel. Through analyses, the order of closing the RCS feedback loops and 
gain changes as necessary will be determined to ensure stability of the model 
during start up of the wind tunnel. Responses will be generated for comparison 
with model responses to monitor the model behavior during the test runs. 


The wind tunnel tests will be accomplished with the sinusoidal gust 
vanes as the model disturbance. Theoretical frequency responses will be generated 
for direct correlation with the wind tunnel results. 
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